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SUMMARY 

The use of a comprehensive analytical model of rotorcraft 
aerodynamics and dynamics is described. This ainalysis is designed 
to calculate rotor performance, loads, and noise; the helicopter 
vibration and gust response; the flight dynamics and handling 
qvialities; and the system aeroelastic stability. The analysis is 
a combination of structural, inertial, and aerodynamic models, that 
ds applicable to a wide range of problems and a wide class of 
vehicles. The analysis is intended for use in the design, testing, 
and evalmtion of rotors and rotorcraft, and to be a basis for 
further development of rotary wing theories. This report describes 
the \ise of the computer program that implements the analysis. 



1 . HK)GRAN SUMMARY 


The computer program calculates the loads and motion of helicopter 
rotors and airframe. First th# trim solution Is obtained; then the flutter, 
flight dynamics, and/or transient behavior can be calculated. Either a 
new job t n be initiated, or further calculations can be performed for an 
old job. 


For a new job, the Input consists of blodc data or an input file 
(the program can create the Input file from the blodc data), and airfoil 
files. Then nuellsts are raeul for additional data, particularly case- 
specific Inputs. One or more cases can be run for a new job. 

For an old job, the input consists of a restart file (written during 
the execution of a previous job), and namelists. Only one case can be run 
for an old job. The job can be resumed either at the point where the trim 
solution was completed, or It can be resumed in one of the subsequent tasks. 
For a trim restart, any or all of the other tasks ca.n be Initiated. For 
flutter, flight dynamics, or transient restarts, only that task can be done. 

For both new and old jobs, a scratch file Is usually needed; and the 
job aay write data on the restart file. In the flutter and flight dynamics 
tasks, eigenvalue data nay be written on a file. 

For both new and old jobs, a case namelist is always read to 
define the job, and a trim namelist is read to define the flight condition 
and analysis tasks. Component and task namelists may be read as required. 

The loads and motion solution is obtained by an iterative process. 

The inner-most loop consists of the rotor and airframe motion calculation, 
for prescribed control positions, Induced velocity distribution, and mean 
shaft motion. Convergence of the motion solution Is determined by comparing 
the calculated harmonics every few revolutions. The next loop consists of 
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the uniform or nonuniform rotor- induced velocity calculation, followed 
hy the motion solution. Convergence is determined by comparing the rotor 
thrust or circulation used to calculate the induced velocity with that 
resulting after the motion has been re-calculated. Before beginning the 
circulation and motion iterations, the blade bending and torsion modes are 
calculated. If the rotor nonuniform induced velocity is used, there is 
an additional cutar loop, consisting of calculation of the rotor wake 
influence coefficients followed by the circulation and motion iterations. 

To calculate the influence coefficients, the prescribed or free wake geometry 
must be evaluated. Having completed the motion solution, the performance, 
loads, vibration, and noise can be evaluated as required. 

The trim analysis proceeds in stages. In the first stage the trim 
solution is obtained for uniform inflow; in the second and third stages the 
trim solution is obtained for nonunlfom inflow, with prescribed or free 
wake geometry respectively. The analysis can stop at any of these stages. 

Within each stage, the aircraft controls and orientation are incremented 
until the equilibrium of forces required for the specified trim state is 
achieved. 

In the flutter analysis, the matrices are constructed that describe 
the linear differential equations of motion, and the equations are analyzed. 
Optionally the equations are reduced to just the aircraft rigid body 
degrees of freedom (by a quasistatic reduction), and the equations are 
analyzed as for the flight dynamics task. 

In the flight dynamics analysis, the stability derivatives are 
calculated and the matrices are constructed that describe the linear differential 
equations of motion. These equations are analyzed (optionally including 
a numerical integration as for the transient analysis). 

In the transient analysis, the rigid body equations of motion are 
numerically integrated, for a prescribed transient gust or control input. 
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2, SUBFROGHAK FUHCTIONS 


The following pages list the subprograms that conSv .ute the 
analysis, and state the priiary function of each subprogram. Only the 
subprograms for rotor #1 are listed; the subprograms for rotor #2 have 
Identical functions. 


Subprogram 

Name 

MAIN Primary job and analysis control 

TIMER Program timer 

IHPTN Inimt for new job 

INPTO Input for old job 

INFTAl Bead airfoil table file 

INPIBI Read rotor xuuiellst 

INPIVI Bead wake namelist 

INFIB Bead body namelist 

INPTLl Bead loads namelist 

ISFTF Read flutter namelist for new job 

INPTS Read flight dynamics naunellst for new job 

INFTT Bead transient namelist for new job 

INPTG Read flutter namelist for old job 

INPTU Re3ui flight dynamics namelist for old job 

INPTV Read transient namelist for old job 

FILEI Read or write Input file 

FUjEJ Read or write trim data file 

FILER Bead or write restart file 

FILEF Read or write flutter restart file 

FILES Read or write flight dynamics restart file 

FILET Reatd or write transient restart file 

FILES Write eigenvalue file 

INIT Initialization 

INITA Initialize environment parameters 

INITC Initialize case parameters 

INITRl Initialize rotor parameters 

INITB Initialize airframe parameters 

INITB Initialize drive train parameters 

CHEKRl Check for fatal errors 
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Subprograa 

Hame 


PRHTJ 

PRNTC 

PRNT 

Himi 

PRITOII 

raNTB 

PRNTP 

PRNTS 

PRNTT 

HINTG 

TRIM 

TRIMI 

raiMP 

FLUT 

FLUTM 

FLUTB 

PLUTRl 

PLUni 

PLUTAl 

FLUTL 

STAB 

STABM 

STABD 

STABB 

STABL 

STABP 

TRAN 

IDAHI 

TRAHP 

TRANC 

CONTRI. 

CUSTU 

GUSTO 

TSBF 

PSRFRl 


Print job input data 

Print ease input data 

Print tria input data 

Print rotor input data 

Print wake input data 

Print body input data 

Print flutter input data 

Print flight dynanlcs input data 

Print transient input data 

Print transient gust and control input data 

Trim 

Calculate trim solution by iteration 
Print tria solution 

Flutter 

Calculate flutter aatrices 

Calculate flutter aircraft aatrices 

Calculate flutter rotor aatrices 

Calculate flutter inertia coefficients 

Calculate flutter aerodynaaic coefficients 

Analyze flutter constant coefficient linear equations 

Flight dynaalcs 

Calculate flight dynaalcs stability derivatives and aatrices 

Print stability derivatives 

Calculate flight dynaalcs equations 

Analyze flight dynaalcs linear equations 

Print flight dynaalcs transient solution 

Transient 

Calculate transient acceleration for nuaerieal integration 

Print transient solution 

Calculate transient gust and control 

Calculate transient control tine history 

Calculate uniform gust tine history 

Calculate convected gust wave shape 

Performance 

Calculate and print rotor performance 
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Subprogru 

Nue 


LOAD 

LOAERl 

LOADHl 

LOADSl 

LOADIl 

LOADP 

LOAIM 

GBOMPl 

FOLRPP 

HISTPP 

NOISRl 

BESSEL 

RAMF 

NODSl 

MOKCl 

HODEBl 

HOlEiG 

HOISAI 

MOIETl 

NOIEKI 

HODBDl 

IHBTCl 

MOEEPl 

Bourc 

BECMC 

MOTHCl 

BOIOMl 

ENGNHl 

WAKEUl 

WAKENl 

INBIMI 

INRTI 

MOTNHl 

MOTNRl 

MOTNBl 

AEROFl 

AEROSl 

AEJROTl 

BODYVl 

EHGHVl 

MOTNFl 

MOTNS 

BODYF 

Bom 


Loads t vibration, and noise 

Calculate and print rotor loads 

Calculate and print hub and control loads 

Calculate and print blade section loads 

Calculate inertia coefficients for section loads 

Calculate fatigue damage 

Calculate mean and half peak-to-peak 

Printer-plot of vake geometry 

Printer- plot of polar plot 

Printer-plot of azimuthal time history 

Calculate and print far field rotational noise 

Calculate J Bessel function 

Calculate rotor/alrframe periodic motion and forces 
Blade nodes 

Initialize blade mode parauaeters 

Calculate blade bending modes 

Calculate Galexkin functions for bending nodes 

Calculate articulated blade flap and lag modes 

Calculate blade torsion modes 

Calculate kinematic pitch-bending coupling 

Calculate blade root geometry 

Calculate blade inertia coefficients 

Print blade modes 

Initialize airframe parameters at trim 
Initialize drive train parameters at trim 
Initialize rotor parameters at trim 
Calculate airframe transfer function matrix 
Calculate drive train transfer function matrix 

Calculate uniform wake- Induced velocity 
Calculate nonuniform wake- induced velocity 

Calculate rotor transfer function matrix 
Calculate inverse of transfer function matrix 
Calculate harmonics of hub motion 

Calculate harmonics of rotor motion 

Calculate blade and hub motion 

Calculate blade aerodynamic forces 

Calculate blade section aerodynamic coefficients 

Interpolate airfoil tables 

Calculate harmonics of airframe motion 
Calculate harmonics of drive train motion 
Calculate rotor generalized forces 
Calculate static elastic motion 
Calculate airframe generalized forces 
Calculate body aerodynamic forces 
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Subprogram 

Name 


VAKECl 

tfAKEBl 

VTXL 

VTXS 

GEONBl 

GEOMRl 

GBOMFl 

MINV 

MINVC 

BIGENJ 

]£RED 

OSITIAN 

GSYSAN 

DElltAN 

SINE 

STATIC 

ZERO 

ZEIBAN 

BOIE 

BOIEPP 

mCKS 

TRCKPP 

GUSTS 

PSYSAN 

DBPRAN 


Calculate Influence coefficients for nonuniform inflow 
Calculate blade insition 

Calculate vortex line segment induced velocity 
Calculate vortex sheet se^ent Induced velocity 

Evaluate wake geometry 

Calculate wake geometry distortion 

Calculate free wake geometry distortion 

Calculate inverse of matrix 

Calculate inverse of complex matrix 

Calculate eigenvalues and eigenvectors of matrix 

Eliminate equations and variables from system of differential 
equations 

Quasistatic reduction of system of linear differential equations 
Analyze system of constant coefficient linear differential 
equations 

Transform equations to state variable form 
Calculate frequency response from matrices 
Calculate static response from matrices 
Calculate zeros 

Transform matrix for zero calculation 
Calculate and printer-plot transfer function (Bode plot) 
Printer-plot transfer function magnitude and phase 
Calculate and printer-plot time history of time- invariant 
system response 
Printer-plot time history 
Calculate and print zms guust response 

Analyze system of periodic coefficient linear differential 
equations 

Transform equations to state variable form 
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3 . NAMELIST, FILE. AND COMMON BLOCK LABELS 


The list below gives the ' ' list labels used by the program, 
and the type of Input data read In each. The corresponding common block 
labels are given in the right-hand column. 


Namelist 


Common Block 

Ubel 


Label 

NLCASE 

Job data 


NLTRIM 

Trim data 

TMDATA 

Numt 

Rotor data 

RIDATA 

NLVAXE 

Wake data 

G1DATA,W1DATA 

NISODf 

Airframe and drive train data 

BDDATA ,BADATA ,ENDATA 

NLLOAD 

Loads data 

LADATA,L1I]ATA 

NLFLUT 

Flutter data 

FLOATA 

NLSTAB 

Plight dynamics data 

STDATA.GCmTA 

NLTRAN 

Tiaasient data 

TNDATA.GCmTA 


The list below gives the files used by the program. The left-hand 


column gives 

the input parameter 

that defines the file unit number. 

Unit 

Number 



NFOAT 

Input data 


NFAFl 

Rotor #1 airfoil 

data 

NFAF2 

Rotor #2 airfoil 

data 

NFRS 

Restart data 


NFBIG 

Eigenvalue data 


NFSCR 

Scratch data 



The list below gives the labels of the comm:nn blocks used by 
the program, and states the type of data contained in each. Only the 
common blocks for rotor #1 are listed; the common blocks for rotor #2 
have identical functions. 
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Conmon Block 
Label 


TMDATA 
RIDATA 
VIDATA 
GIDATA 
BDDATA 
BADATA 
SNDATA 
LI DATA 
UDATA 
GCDATA 
TNLATA 
STDATA 
FLDATA 
AITABL 
UNITNO 
CASECM 

TRIMCM 

RTOIGM 

RHICM 

BODYCM 

BNGNCM 

GUSTCM 

CONTCM 

CONVCM 

MDICM 

INCICM 

WKVICM 

MNHICM 

ABSICM 

MNRICM 

MNSCM 

AEPICM 

QRICM 

QBDCM 

WGICM 

WKClCM 

ASMNCN 

LDMNCM 

FLMCM 

FLMICM 

FLMACM 

FLINCM 

FLAECM 

STDCM 

STMCM 

TRAHCM 


Input trim data 

Input rotor data 

Input wake data 

Input free wake geonetry data 

Input airframe data 

Input airframe aerodynamics data 

Input drive train data 

Input rotor loads data 

Input airframe loads data 

Input gust and control data 

Input transient data 

Input flight dynamics data 

Input flutter data 

Rotor airfoil tables 

Input/outpuc unit numbers 

Job description 

Calculated trim data 
Calculated rotor data 
Transfer fvinction matrices 
Calculated airframe data 
Calcvilated drive train data 
Gust and transient control 
Aircraft controls and motion 
Circulation and motion convergence 
Blade modes 

Rotor inertial coefficients 
Induced velocity 
Hub motion 

Blade section aerodynamics 
Rotor motion and airframe vibration 
Static elastic motion 
Rotor forces 

Rotor generalized forces 
Airframe generalized forces 
Hake geometry 

Wake influence coefficients 

Calculated motion data 
Calculated loads data 
Flutter matrices 
Flutter rotor matrices 
Flutter airframe matrices 
Flutter inertial coefficients 
Flutter aerodynamic coefficients 
Flight dynamics stability derivatives 
Flight dynamics matrices 
Calculated transient data 
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4. PROGRAM SKELETON 


The following pages present a schematic of the program, showing 
the haslc flow of control and the major loops, options, and decisions. 

The appearance of a subprogram name (always In capital letters) means 
that the subprogram Is called at that point In the analysis. The contents 
of a subprogram are shown by means of a thiree-slded box appearing below 
the subprogram name. The schematic defines the Input and output structure 
of the program. Timer calls and trace-debug prints are also shown. 
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KAIN 


read naunellst NLCASB 
if new job and BLKDAT > 0 
DATS (for FILEID} 

TIME (for FILEID) 

FILEI (input file write) 

PRNTJ 

for JCASE - 1 to NCASBS 

TIMER (initialize) 

TIMER 

DATE (for IIENT) 

TIME (for IDENT) 
if new job 

INPTN 

INIT 


INITA 

INITC 

INITRl 

INITR2 

INITB 

INITE 

CHBKRl 

CHEKR2 


if old job 

INPTO 

PRNTC 

if new job or trim ret, cart 
TRIM 

FILEJ (trim data scratch file write) 
ANTYPE(l) 0 or flutter restart 
FLUT 

FILEJ (trim data scratch file read) 
if ANTYPE(2) f 0 ot flight dynamics restart 
STAB 

FILEJ (trim data scratch file read) 
if ANTTFE(3) / 0 or transient restart 
TRAN 

TIMER 

TIMER (print) 
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IHPTH 

PILBI (input file read) BLKD/^T.HUFILE 

read naaelist NLTRIN 
if OPRBAD(l) ^ 0 

INPTRl 

read naaelist HlflTR 

if 0FREAD(2) / 0 

INPTWl 

read naaelist NLWAKS 

if 0HIBAD(3) ^ 0 
IKPTB2 

|read naaelist NLRTO 

if 0FREAD(4) / 0 
IHPTO2 

|read naaelist HLHAKE 

if 0PREAD(5) / 0 
IHPTB 

jread naaelist NLBODY 

if 0PREAD(6) * 0 
INPTLl 

[read naaelist NUiOAD 

if OPREAD(7) I* 0 

INPTL2 

[read naaelist NLLOAD 

if OPREAD(0) a 0 
IKPTF 

read namelist NLPLUT 

if 0FRBAD(9) f 0 
INPTS 

jreaui namelist HLSTAB 

if 0PREAD(10) I* 0 
INPTT 

read namelist NLIBAN 

if first case 

INPTAl 

[read airfoil file 
INPTA2 

(read airfoil #2 file 
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TRIM 

TIMER 

if trim restart, go to restart entry point 

unifom inflow 
if ITERU / 0 

TRIMI 

if HPRNTT - 1 
lERF 
LOAD 

nonunifom inflow, prescribed wake geometry 
for IT = 1 to ITERR 
WAKECl 
VAKEC2 
TRIMI 

if IT = multiple NPRNTT 
PERF 
LOAD 

nonunifom inflow, free wake geometry 
for IT = 1 to ITKRF 
WAKECl 
WAKEC2 
TRIMI 

if IT = multiple NPRNTT 
PERF 
LOAD 

FILER (restart file write) 

Ifilei 

FILEJ 

trim restart entry point 
PRNT 

PRNTC 

if NPRNTI i 0 
PRNTRl 
PRNWl 
PRNTR2 
PRNW2 
IRNTB 

MOEBPl 

MOEEP2 

TRIMP 

lERF 

LOAD 

TIMER 


NPRNTP 

NPRNTL 


LEVEL(l) 

LEVEL(2) 


NPRNTP 

NPRNTL 


LEVEL(l) 

LEVEL(2) 


NPRNTP > 0 
NPRNTL > 0 

RSWRT / 0 
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A A A\ A A AV Av 



ORtGtHAi PACE 



’SSii 9 

RAMF 

if MTOIM 0 or no tri* iteration, return 
if DEBUG(4)» 1, print tri« iteration 
for COUNTT = 1 to MTRIM 

if COUNTT-1 » multiple MTOIMD, construct D"! 
for I = 1 to MT 

increment controls 
RAMP 

MIHV 


increment controls 
HAHF 

if BEBUG(4)> 1, print trim iteration 
test trim convergence 


OPTRIK 

CPTOIM 

KPTRIM.CPTRIM 


i«RF 


TIMER 

iSRFRl 

IVRFR2 

TIMER 
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LO^D 


(tjMER 
]l Ami 
\0Am2 
I »*IHSR 


LOAmi 


1 )TRB1 

« ’ MALOAD i 0 
GBOMEl 
HISTPP 
GEOMPl 
POUiPP 
HISTPP 

If HHLOAD ^ 0 
LOAMI 


LOADf 

LOAIF 

HISTPP 


for IB - 1 to MRLOAD 
LDAmi 


LOADTl 
LO* H 
LOAIF 
HISTPP 


for IN * 1 to MHDISE 
HolSNl 



NPLOT(l-4) 

NUKGMP 

NPL0T(5-67) 

NPL0T(5-67) 


NPLOT( 68-71) 


NPL0T(72-75) 


- 16 - 



FLUT 

TIMra 

for OPFLOW ^ 0 (constant coefficients) 

if flutter restart, go to restart entry point 
FLUTM 

FILEF (restart file write) 

flutter restart entry point 

PRNTF 

MODEPl 

MOOBP2 

FLUTL 


TIMER 

CSYSAM 

PILES (elg-jnvalue file write) 

BOSE 

TRACKS 

GUSTS 

TIMER 

if OFFDAM f 0 
STABD 
STABS 

for OPFLOW> 0 (periodic coefficients) 
for NT = 0 to MPSIPC 
FLUTM 

if NT - MPSIPC 
IRNTF 
MODEPl 
N0DBP2 

PSYSAN 

If NT - MPSIPC 

FILES (eigenvalue file write) 


RSWRT / 0 


ANTYPE(l) f 0 

ANTYPE(2} / 0 
ANTnE(3) I* 0 
ANTYPE(4) / 0 
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FLUTM 


KOIXl 

M0DS2 

FLUTSl 

FLUTR2 

FLUTB 


BODYF 


FLUTRl 

NB - NBLABE If OPFLOW >0, I if OPFLOM - 0, MPSIGC if OPFLOW < 0 
for JPSI - 1 to NB 
FLUTIl 

FLUTAl 

Ifor IR - 1 to MRA 
AEROSl 
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STAB 


TIMKR 

PRHTS 

If flight dynamics restart, go to restaort entry point 
STABM 


for ID - 1 to 21 

increment controls or notion 
for IT ” 1 to ITSRS 


WAKBCl LEVEL(1^ 
WAKBC2 LEVEL(2) 
RAMF 

PBRF KPRNTP 
LOAD NHJNTL 


FILES (restart file write) RSWRT i 0 


flight dynamics restart entry point 

STABD 

STABS 

TIMER 


STABE 


1 for IRQ - 1 to 12 
IXRED 
STABL 


TIMER 

CSYSAH 

FILES (eigenvalue file write) 

IBODB 

TRACKS 

GUSTS 

numerical integration 
MINV 
STABP 


PRNTG 


for IT 


1 to tmax/tstep 

TRANC 


OONTRL 

IGUSTU 

GUSTC 


if IT - multiple NIRHTT 
STABP 


TRCKPP 


TIMER 


EQTYPE(IEQ) f 0 


ANTYPEd) / 0 

ANTYPE(2) ¥ 0 
ANTYPEp) y 0 
antype( 4) y 0 
ANTYPE(5) y 0 


OPTRAN » 1 
OPTRAN - 2 
OPTRAN - 3 
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vv nVv 



TRAN 


TIMER 

FRMTT 

|frntg 

If txanslent restart, go to restart entry point 

MINV 

TRAMP 

for IT - 1 to TMAX/tSTBP 
TRANG 


CONTRL 

GUSTU 

GUSTO 


TRANI 

for IT “ 1 to ITRRT 
VAKECl 
WAKEC2 
RANF 


If IT - Bultiple NPRHTT 
TRANP 
PERF 
LOAD 

if IT - Multiple NRSTRT 

FILET (restart file tirite) 
transient rest^t entry point 

TRCKFP 

TIMER 


OPTRAN « 1 
OPTRAN = 2 
OPTRAN - 3 


LBVBL(l) ^ 1 
LBVEL(2) 1 


NPRNTP> 0 
NFRNTL> 0 

RSWRT / 0 
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RAMF 

TIMER 

BODYC 

f«)TNCl 

MOIEl 

BODYMl 

IlNRTI 


M0TNG2 

M0DB2 

BODYM2 

jiMRTI 

for OOUNTC = 1 to ITERC (circulation iteration) 
WAKEUl 
WAKENl 
WAKEU2 
WAKEN2 

for OOUNTM =» 1 to ITERM (motion iteration) 
INRTMl 


INRTI 


INRTM2 


INRTI 


BNGNC 

ENGNMl 

IinrtT 

BNGNM2 

I INRTI 

for JPSI - 0 to MREV*MPSI by MPSIR (V loop) 
MOTNHl 
MOTHRl 
M0TNH2 
M0TNR2 
BODYVl 
ENGNVl 
MOTNFl 
B0DYV2 
ENGNV2 
M0TNF2 
MOTNS 

test notion convergence 
test circulation convergence 

BODYF 

|bodya 

TIMER 


EPMOTN 

EPGIRC 
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MOISI 


TIMER 

MOIECl 

if A» > EFNOIE 

MOEEBl HINGE f 2 

Imoibg 

MINV 

EIGENJ 

HODEAl HIIKJE = 2 

MOHK1 

HOIXDI 

NOIETI 


MINV 

EIGENJ 

INRTCl 

TIMER 


MOTNRl 

TIMER 

for JP “ JPSI + 1 to JPSI + MPSIR (Y step) 
MOTNBl 
AEROFl 


Ifor IR * 

1 to MRA 

1 AEROSl 


AEROTl 


TIMER 
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origin 
OF POOR 


page is 

Q\3AEn:Y. 


WAK5C1 



LEVEL = 2 

DBV> 0. 

IBV:» 0. 

IN'. L0W(3) * 3 
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CSYSAH 


DmAN 

EIGENJ 

SINE 

STATIC 

ZERO 


ZETRAN 

EIGENJ 


BODE 


1£TRAN 

EIGENJ 

ZERO 


IZETRAN 

{EIGEN J 

BO DEPP*" 


mCKS 


DEIBAN 

EIGENJ 

MINVC 

TRGKPP 


GUSTS 


DETRAN 

EIGENJ 

MINVC 


PSYSAN 


DBFRAN 

EIGENJ 



5 . JOB STRUC1UHE 


In this section the structure of a job to run the program is 
defined. The basic structure consiets of the folloHing ste-s: 
l) File definition as required for job 
?.) Block data load for airframe and each rotor 

3) Main program call 

4) Namelist 4NLCASE 

5) Namelist &NLTRIM (for each case) 

6) Component and task namelists as required 
File definition parameters j 

a) RET = T Erase file at logoff 

b) DISP = NEW New file to be created 

c) DISP = OLD Existing file 

Sample jobs are presented below. 


New job, 2 cases; trim analysis; block data input, basic namelist input, 
same airfoil table for both rotors 

DDE? FT50F001 , , SCRATCH, DISP=*NEW,RE'I^T 
DDBF FT41F001,, AIRFOIL, DISP=OLD 
LOAD HELA; LOAD HELRl; LOAD HE1R2 
CALL MAINPROG 

&NLCASE JOB-0 , NCA2ES»2,RSWRT-0,BLKDAT»-1 , 

NFAFl-41 ,NFAF2-41 ,NFSCR-50,NFRS--1 ,NFEIG»'-1 , 

&IND 

ANLIRIM VKTS=x, ,v/DLL“X. ,LATCYOx. ,LNGCYC«x. ,PBnAL*x. ,APITCH»x. ,AROLL-x. , 
ANTyra-3*0 ,0I®EAD-10*0 , 

&SND 

&NLTRIM data for second case,4END 
^BND 
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New job, 1 case; trin, flutter, flight dynamics, and transient analysis; 
blodc data input, all namelist inputs, different airfoil table for each 
rotor; write eigenvalue file 

DEBF FT50F001,, SCRATCH, DISP-HBf,RrP-T 
DIBF miFOOl , , AIRPOILl , DISIHDLD 
DEBF FTbZFOOl , .AIRF0IL2, DISIH3LD 
DI*F FT45F001 , ,BIGBK,DISP-HEN 
LOAD HBLA; LOAD HBIiil; LOAD 
CALL MAINFROG 

mCASE JOB-0,NCASES>l,RSVRT^0,BUCDfl'iy-l, 

HFAFl-41 ,NPAF2-42,NFSCR«50,NFRS--1 ,NFBIG*45, 

AEND 

&HLT xM VKTS-x., 

G0LL«x. ,LATCYOx. ,LNGCYOx. ,PBDAL-x, ,APITCH-x. ,AR0LL=x. , 
ANTyPE-3*l ,0HIEAI>«10*1 , 

ABND 

AHIitlR data,AEND 
ANLHAKE data,AEND 
AMLRTR data,AEND 
ANLWAKE data,AEND 
ANLBODY data,AEMD 
ANLLOAD data,&rAD 
ANLLOAD data,AEMD 
&NLFLUT data,ASND 
ANLSTAB data,AEND 
ANLTRAN data &END 
^ND 


New job, 1 case: trim analysis; block data input and write input file 

DDEF FT50F001,,SCRATCH,DISP“NEM,RBT-T 
DDBF FT41F001 , , AIRFOIL , DISP-OLD 
deep FT40F001 , , INPUT, DISP-NBH 
LOAD HELA; LOAD HBLRl: LOAD HELR2 
CALL MAINITOG 

&NLCASE JOB-O , NCASB-1 ,RSWR‘I^0 .BLKnAI^l , 

NFAFl-41 , NFAF2-41 , NFSCR-50 , NPRS=>- 1 , NFEIG--1 , NFDA'1^40 , 

AEND 

ANLTRIM data,AEKD 
^ND 
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New job, 1 case; trial analysis] read input file 

DIEF FT50F001 , , SCRATCH , DISP-NBN ,RET>»T 
DIEF FT41F001,,AIRF0IL,DISIM)LD 
DDEF FT40F001,, INPUT, DISPK)LD 
CALL MAINHiOG 

&NLCASE JOB-0 , NCASES«1 ,RSWRT“0 ,BLKDAT“0 ,RIFILB=1 , 
NFAFl-41 ,NFAF2»41 ,NFSCR-50,NFRS=-1 ,NFBIG*-1 ,NFDAT“40 , 
&END 

ANLliaM data,&END 
^ND 


New job, 2 cases; trin and flutter analysis; write restart file 

DIBF FT50F001 , , SCRATCH , DISP=NBH ,RET»T 
EDBF FT41F001,, AIRFOIL, DISIK)LD 
DIEF FT44F001 , .RESTARTl , DISP=NBif 
DDEF FT44F002,,RESTART2,DISP«NEif 
LOAD HELA; LOAD HEUa; LOAD HEIB2 
CALL MAINIROG 

&NLCASB J0B=0,NCASBS=2,RSWR‘I>=l,BLKDA'iy-l, 

NFAF1=41 , NFAF2=41 , NFSCR=50 , NFEIG*-1 , NFRS»44 , 

AEND 

&NLTRIM data for first case, 

ANTYPE=1 , 0 , 0 ,0PREAD{ 8)-1 , 

AEND 

ANLFLUT data,AEND 
ANL'TOIM data for second case.AEND 
ANLFLUT data,ASND 
%BND 


Old job; trin restart with flutter aralysis 

DIEF FT50F001 , ,SCRATCH, DISI^HEH ,RET^T 
DDEF FT44F001,, RESTART, DISP-OLD 
CALL NAINPROG 

ANLCASE JOB-1 ,RSWRT-1 , START-1 , 
NFSCR-50,NFEIG--1 ,NFRS-44, 

AEND 

ANLTRIM ANTYPE-1,0,0,OPREAD(8)-1, 
AEND 

ANLFLUT data, AEND 
;;aEND 
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Old job; flutter restart 

DEBF FT50F001 , .SCRATCH, DIS>MW ,RTM 
DDBF PT‘|4F001,,R1START,IIIS1M)LD 
CALL HAINHiOG 

ANLCASB JOB'l.RSHRT-O.STARl^a. 
MFSCR-50,!fFSIG--l .HFRS-44, 

ABMO 

AMLTRIN OHtEAD(8)-l, 

ABMD 

AMIfLUT data.AEMD 
^CBIID 
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6. INPUT EESCRIPTICH 


In this section the input variables for the progran are defined. 
The variables are catagorized according to the nanelist that reads them. 
The program namelist labels are listed in the table below. 


Namelist 

Label 

NLCASS Job data 

NLTRIM Trim data 

NLRTR Rotor data 

NLHAKE Hake data 

N1£0DI Airframe and drive train data 

NLLOAD Loads data 

NIFLUT Flutter data 

MLSTAB Flight dynamics data 

NLTRAN Transient data 


The corresponding common block labels, for the blodc data form of input, 
may be obtained from Section 3* In the description of the input parameters 
for the rotor, the variables NBH and Nllf are usedt 

a) NBN is the index of the highest-frequency blade bending 
mode used in the analysis; 

b) NTM is the index of the highest-frequency blade torsion 
node used in the analysis. 
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M»«ellst HLCASE 


JCB 

RSVRT 


NCASES 

BLKD^I 


RIFILE 


START 


HFIAT 

NFAFl 

NFAF2 

NFRS 

NFEIG 

NFSCR 

NUIN 

NUOUT 

NUIB 

NUPP 

NULIN 


Integer paiaeeter defining jobt BQ 0 for new job (default); 
ME 0 for old job or restart (one case only) 

integer paraaeter controlling restart file write; 0 to 
suppress write (defai0.t) 

Mew job only 

number of cases (default ~ l) 

Integer paraaeter defining input source: 

B(1 0 reaul input file (default) 

GT 0 use loaded blodc data and write input file 

LT 0 use loaded blodc data 

integer paraaeter controlling input file read: 

E(1 0 read file for first case only 

HE 0 read file for every case (default) 

Cld job only 

integer paraaeter defining task: 

1 for tria restart (default) 

2 for flutter restart 

3 for flight dynamics restart 

4 for transient restart 

tria restart can be followed by any or all of the other tasks 
(as defined by AMTYPE); for flutter » flight dynamics, or 
transient restart, only that task can be done 


Input/output unit nuabers 

input data file (new job only); defaiQt = 40 

rotor #1 airfoil file (new job only); default * 41 

rotor #2 airfoil file (new job only; only if have two rotors); 
default ■= 42 

restart file (no file write if LE O); defa'ilt * 44 

eigenvalue file (no file write if LE O); default “ 45 

scratch file; default =■ 50 

namelist input; default 5 

printer (and debug level l); default * 6 

debug output (levels 2 and 3) 5 default * 6 

printer- plots; default ■ 6 

linear system analysis; default « 6 
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Naaellst NLTOIM 


OraEAD(lO) integer vector defining namelist read structure; SQ 0 

to suppress read: 

components (new job only) 

(1) NLRTR, rotor #1 

(2) NLWAKE, rotor 

(3) NLRTR, rotor UZ 

(4) NLWAKE, rotor #2 

(5) NLBODY 

tasks 

(6) NLLCAD, rotor 

(7) NLLOAD, rotor #2 

(8) NLFLUT 

(9) NLSTAB 

(10) NLTRAN 

NPRNTI integer parameter controlling input data print; £Q 0 

for short form print only 

ANTYPE(3) integer vector defining tasks for new job or trim 

restart; EXi 0 to suppress: 

( 1 ^ flutter 

(2; flight dynamics 

(3) transient 

TITLE(20) title for job and case (80 characters) 

COLE alphanumeric code for job and case identification; 

4 characters 

OPUNIT integer parameter designating unit system: 1 for 

English units (ft-slug-sec) ; 2 for metric units (m-kg-sec) 

NROTOR number of rotors 
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NLTRIM 


DBBUG(25) 


Integer vector controlling debug prints 

0 no debug print 

1 trace print 

2 low level print 

3 high level print 


( 1 ) 

(2J 

(3 

\i 

( 6 ) 

(7 

8 

(9! 


(ii; 

(12 

(13) 



tine (sec) at which debug print enabled 
input, 2-3 (iNPTx) 

initialization, 2 (INITC, INITR, INITB, INITE) 

trim iteration, 1-2 (TRIMi) 

loads, 2 (LOADI) 

flutter matrices, 2-3 (FLUTW) 

flutter coefficients, 2-3 (FLUTI, FLUTA) 

flight dynamics, 2-3 (STABM, STABE) 

transient, 2 ( IRANI ) 

rotor/airframe motion and forces, 2-3 (RAMF) 

blade modes, 2 (NOIS, MO^x) 

inertia coefficients, 2 (INRTC) 

airfr 2 une constants and matrices, 2 (BODYC, ENGNC, 

MOTNC, BO Dim, ENGNM) 

induced velocity, 2 (WAKEU, WAKEN) 

rotor matrices, 2-3 (INRW) 

hub/airframe motion and generalized forces, 2 

(MOTNH, BODYV, BNGHV, MOIRF, MOTNS) 

rotor motion, 2-3 (MOTNR) 

rotor aerodynamics, 2-3 (AEROF) 

blade section aerodynamics, 3 (AERCS) 

body forces and aerodynamics, 2 (BODYF) 

wake influence coefficients, 2 (WAKEC) 

vortex line and sheet, 3 (VTXL, VTXS^ 

prescribed wake geometry, 2-3 (GBOMR) 

free wake geometry, 1-3 (GEOMF) 

timer, 1 (TIMER) 
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VKTS 

VSL 

'■TIP 

RPM 


CPDENS 


ALWSL 

TEMP 

DENSE 

OPGRND 

HAGL 

OFENGN 


AFLAP 

RTURN 


original page is 
OE POOR QUAUTXi 

NLTRIM 

aircraft speed Y (knots) 
velocity ratio V/^R 

input either VBL or VKTS by namelist; if neither 
parameter is defined, V = 0 is used 

rotor tip speed .T2.R (ft/sec or m/sec) 

rotor #1 rotational speed (rpm) 

input either VTIP or RPM by namelist; if neither 
parameter is defined, the normal tip speed VTIPN 
is used; rotor #2 speed is calculated from the 
gear ratio T3ATI0 

integer parameter defining specification of aerodyptriic 
environment! if 1, given altitude and standard day; 
if 2, given altitude and temperature; if 3, given density 
and temperature 

altitude above mean sea level (ft or m), for CPDENS = i or 2 

air temperature (°F or °C), for OPDEIG = 2 or 3 

air density (slug/ft^ or kg/m^), for OPIENS = 3 

integer parcimeter controlling ground effect analysis; 

EQ 0 for out of ground effect, NE 0 for In ground effect 

altitude helicopter center of gravity above ground for 
ground effect analysis (ft or m) 

integer parameter specifying engine state: 1 for autorotation 

(engine inertia, engine damping, and throttle control torque 
zero: no engine speed degree of freedom); 2 for engine out 
(engine damping and throttle control torque zero); 0 for 
normal operation 

wing flap angle Sp (deg) 

for free flight, trim turn rate Vp (deg/sec), positive to right 
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NLTRIM 


COLL 

LATCYC 

LNGCYC 

PEDAL 

APITCH 

AROLL 


ACLIMB 

AYAW 


KPSI 

MPSIR 

MREV 

ITERM 

SPMOTN 

ITERC 

EPCIRC 


initial values of controls (trimmed as appropriate) 

collective stick displacement £ or AQ (deg), 
positive up ° govr 

lateral cyclic stidc displacement J (deg), positive left 

longitudinal cyclic stick displacement C (deg), positive 
aft ® 


pedal displacement (deg), positive to right 

for free fll^t, aircraft pitch angle (deg), positive 

nose up: for wind tunnel, rotor shaft angle of attack 
(deg), positive nose up 


for free flight, aircraft roll angle (deg), positive 

to right 


( ©p and define orientation of body axes relative to 

earw axes) 


for free flight, aircraft climb angle ©pp (deg), positive up 

for free flight, aircraft yaw angle H'pp (deg), positive to 
right: for wind tunnel, test module yaw angle 
positive to right 


(©pp and Vpp define orientation of velocity axes relative 
to eSrth axes: = V sin ©pp and = V slnVpp cos©pp) 


number of azimuth steps per revolution in motion and loads 
analysis, maximum 36 : for nonuniform inflow must be multiple 
of number of blades: for free Wcike geometry, maximum 24 

in harmonic motion solution, number of azimuth steps between 
update of airframe vibration and rotor matrices 

in harmonic motion solution, number of revolutions between 
tests for motion convergence 

maximum number of motion iterations 

toleirance for motion convert :ence (deg) 

maximum n\imber of circulation iterations 

tolerance for circulation convergence ( ) 
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D0F(54) 


D0FT(8) 


MHARM(2) 


MHARMF(2) 


LEVEL(2) 


ORIGINAL’ PAGE IS 
OF POOR QUALITY 

NLTRIM 

integer vector defining degrees of freedom used in 


vibratory motion solution, 0 

if not used; order: 


rotor #1 

^1 ••• ^10 

Pq ••• P4 

Pc 

rotor #2 

0 

m 

• 

Pq ••• P4 

Pc 


(bending, max 10 ) 

(torsion, max 5) 

(gimbal/ teeter) 


airframe Xp Zp 

(rigid body) 

drive train *P 

S J. 6 

(rotor/engine speed) 


•Is'p • • • 

(flexible body, max lO) 

t govri govT2 
( governor) 


integer vector defining blade bending degrees of freedom 
used for mean deflection (subset of DOF), 0 if not used; 
order t 

rotor #1 qg <14 

rotor #2 q^ qg q^ q^. 

(bending, max 4) 


number of hamonics in rotor motion analysis; maximum 20; 
EQ 0 for mean only 

(1) rotor #1 

(2) rotor #2 

number of harmonics in airframe vibration analysis 
(harmonics of N/rev); maximum 10; EQ 0 for static elastic 
only; suggest LE MHARM/NBLADE, and the same value for 
both rotors if coupled hub vibration used (see 0PH\^IB) 

(1) rotor #1 

(2) rotor #2 

integer parameter specifying rotor wake analysis level: 

0 for uniform inflow, 1 for nonuniform Inflow with 
prescribed wake geometry, 2 for nonuniform Inflow with 
free wake geometry (must be consistent with INFLOW) 

(1) rotor #1 

(2) rotor #2 
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NiraiM 


ITBRU 

IlBRR 

ITERF 

NHINTT 

NPRNTP 

NHUnii 

MTMM 

MTRIMD 

DELTA 

FACTOR 

EPTRIM 

OPGOVT 


CXTRIM 

xraiM 

CTTRIM 

CPTRIM 

CYTRIM 

BCTRIM 

BSTRIM 


nuaber of wake and trim iterations 

at vtnifoxm inflow level; BQ 0 to skip 

at nonuniform inflow/prescribed wake geomet-ry level; 

EQ 0 to skip 

at nonuniform inflow/free wake geometry level 

integer parameter ni trim/performance/load print 
every n-th iteration; LS 0 to suppress 

integer parameter controlling performance print; LE 0 to 
suppress 

integer parameter controlling loads print; LE 0 to suppress 

maximum number of iterations on controls to achieve trim 

number of trim iterations between update of trim derivative 
matrix 

control step in trim derivative calculation (stick displacement t 
deg) 

factor reducing control Increment in order to improve trim 
convergence (typically 0.5) 

tolerance on trim convergence 

integer parameter specifying governor trim 

0 trim collective stick Sg 

1 trim rotor #1 governor 

2 trim rotor #2 governor 

3 trim both rotor governors 

targets for wind tunnel trim cases 

X/q (ft^ or m^) 
or cyv- 

Cp/r 

Cy/O* 

(deg) 

(deg) 
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OPTRIM 


NLTRIM 


Integer pareuneter specifying trim option 
free flight cases 


OPIRIM = 0 

no trim 


1 

trim 

forces and moments with 



«^o 

Sc Sg Sp 


2 

trim 

forces and moments with 



So 

Sc Sg &p H'pp 


3 

trim 

forces, moments, and power 

with 


So 

Sc Ss Sp ^FT ^FT ^FP 



trim 

forces, moments, and power 

with 


So 

Sc Sg Sp ®FT **^FP ^F? 



trim 

symmetric forces and moments with 


So 

Ss ^FT 


6 

trim 

symmetric forces, moments, 

and power 


So 

Sg ®ft ^fp 



wind tunnel 

cases 








OPTRIM = 10 

no trim 







11 

trim 



With 





12 

trim 

V- 


with 

®T 




13 

trim 

Gp/e- 


with 

^0 




14 

trim 

Ps 

with 

K 

Ss 



15 

trim 

Cp/v 

Pc Ps 

with 

So 

Sc 

Ss 


16 

trim 

Cj^/y 

CpA- 

with 

So 

Sc 

Ss 


17 

trim 

Ci/y 

C/c- Cp/*- 

with 

So 

Sc 

©T 


18 

trim 


Pc Ps 

with 

So 

Sc 

Ss 


19 

trim 


X/q Cy/ct- 

with 

So 

Sc 

Ss 


20 

trim 

Cj/y 

X/q C^/*- 

with 

So 

Sc 

Or 


21 

trim 

V"- 

'‘/l fc Ps 

with 

^0 

Sc 

Ss 

«T 

22 

trim 

^c 


with 

^6 




23 

trim 

Cp/«- 

p= 

with 

So 

Ss 



24 

trim 


0/<r 

with 

So 

Ss 



25 

trim 


y<r 

with 

So 

St 



26 

trim 


V'" P= 

v:ith 

So 

Ss 



27 

trim 

Ci/c- 

X/q 

with 

So 

Ss 



28 

trim 

yr 

X/q 

with 

Sb 




29 

trim 


x/q (J, 

with 

So 

Ss 

®T 
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NT.T?IM 


WEIGHT see namelist HLBODY 
IXX 

m 

IZZ 

IXI 

IXZ 

lYZ 

ATILT 

PSCG 

BLCG 

WLCG ' ' 
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original page is 

QP POOR 


Namelist NLRTR 


TITLE(20) 

TYPE 

VTIPN 

RADIUS 

SIGMA 

GAMMA 


NBLAIE 

TDAMPO 

TDAMPC 

TDAMPR 

NUGC 

NUGS 

GDAMPC 

GDAMPS 

LDAMPC 


LDAMPM 

LDAMPR 

GSB(NBM) 

GST(NIM) 

ROTATE 


title for rotor and waJce data (80 characters) 

rotor Identification (4 characters); suggest MAIN, FRNT, 
or RGHT for rotor #1; and TAIL, HEAR, or LEFT for rotor #2 

normal tip speed ^ZR^ (ft/sec or m/sec) 

olade radius R (ft or m) 

solidity ratio ^ = Nc /rtR (based on mean chord) 

blade Lock n’jmber ^ ac R /l-^ (based on standard 

density, a « 5*7» anS mean°clwrd) 

( ^ and r are only used *o calculate the normalization 
parameters c and I. ) 

m D 

number of blades 

control system damping (ft-lb/rad/sec or m-N/rad/sec) 
collective 

cyclic 

rotating 

longitudinal glmbal natural frequency or teeter natural 
frequency (!>er rev at normal tip speed VTIPN) 

lateral glmbal natural frequency (per rev at normal tip 

speed VTIPN) 

longitudinal glmbal daunplng or teeter damping C„ 
(ft-lb/rad/sec or m-N/rad/sec; ^ 

lateral gi. bal damping C^g (ft-lb/rad/sec or m-N/rad/sec) 

linear lag damper coefficient (ft-lb/rad/sec or 
m-N/rad/sec); estimated damping if a nonlinear damper is 
used f LDAMPM GT 0.); the lag mode has structural damping 
also (GSB) 

maximum moment of nonlinear lag damper, M.^^ (ft- lb or m-N); 
linear lag dami>er used If LDAMPM EQ 0. 

lag velocity 5 tq where maximum moment of lag damper occurs 
(rad/sec); hydraulic damping below ^■nd friction damping 

above 

bending mode structural damping g_ 
tcrcicn mode structural damping g^ 

Integer parameter specifying rotor rotation direction: 

1 for counter-clockwise, -1 for clockwise (viewed from 
above) 
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NLRTR 


0PHVIB(3) 

BTIP 

OPTIP 

LINTO 

TWISTL 

OPUSLD 

OPCX)MP 

INFL0W(6) 


RROCT 

RGMAX 


integer parameter controlling hub vibration contributions; 
gravity and static velocity terms always retained; 0 to 
suppress t 

(ll vibration due to this rotor 

(2) vibration due to other rotor (must 

supress if ^ 

(3) static elastic motion 

tip loss parameter B 

integer parameter specifying tip loss type: 1 for tip loss 

factor, 2 for Prandtl function 

integer parameter specifying twist type: EQ 0 for nonlinear 

twist, NE 0 for linear twist 

linear twist rate (deg); used to calculate TWISTA and 

WISTI if LINTW NE 0™ 

integer parameter controlling use of unsteady lift, moment, 
and circulation terms: if 0, suppress; if 1; include; if 2, 

zero for stall ( 15 ® <. |b<\ < 165°) 

integer parauneter controlling aerodynamic model, EIQ 0 for 
incompressible loads 


Inflow model 

integer vector defining Induced velocity calculation 
(must be consistent with LEVEL) 

(1) at this rotor: 0 for uniform, 1 for nonuniform 

( 2 ) at other rotor: 0 for zero, 1 for empirical, 

2 for average at hub, 3 for nonuniform (only 

if = 1 ) 

( 3 ) at wing-body: 0 for zero, 1 for empirical, 

2 for nonuniform 

( 4 ) at orizontal tail: 0 for zero, 1 for 

empirical, 2 for nonuniform 

( 5 ) at vertical tail: 0 for zero, 1 for empirical, 

2 for nonuniform 

(6) at point off rotor disk: 0 for zero, 1 for 

nonuniform 

root vortex position for wake model, 

tq /r (induced velocity calculated using maximum bound 
circulrtion magnitude outboard of 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


NLRTR 


NRA 

RAB(HRA+1) 


CH0RD(MRA) 

XA(MRA) 

thetzl(mra) 


TVISTA(MRA) 

xac(mra) 

HOORRL(MRA) 

MCORRD(MRA) 

MCORRM(MRA) 


Blade section aerodynanlc characteristics 

nunber of aerodynaaic segments; maxiauni 30 

radial stations r^ at edges of aerodynanlc segments; 
sequential, from root to tip 

Following quantities are specified at midpoint of 
aerodynaunic segment 

blade chord, c/R 

offset of aerodynamic center aft of elastic axis, x^/R; 
X. is the point about which the moment data in the 
airfoil tables is given 

incremental pitch of aero lift line, can be 

included in TWISTA; Pitch of the axis 

corresponding to zero angle of attack in the airfoil 
tables, relative to the twist angle (WISTA) 

blade twist relative .75Rt (deg) 

offset of aerodynamic center (for luustcady aerodynamics) 
aft of elastic axis, 

Mach number correction factor lift 

Mach nunber correction factor fj^ « drag 

Mach nunber correction factor f^ “ M 


efr ' 


M for moment 


MRI 

RI(MRI) 

MASS(MRI) 

EIXX(MRI) 

EIZZ(MRI) 

XI(MRI) 

XC(MRI) 

KP2(MRI) 

ITHETA(MRJ) 

GJ(MRI) 

IVISTI(MRI) 


Blade section inertial and structuial characteristics 

number of radial stations where characteristics defined; 
maucimum 31 

radial stations r/R; sequential, from root to tip, 

RI(1) = 0. and Rl(MRl) •= 1. 

section mass, m (slug/ft or kg/m) 

chordwise bending stiffness (ib-ft or N-m ) 

flapwise bending stiffness (ib-ft^ or N-m^) 

offset of center of gravity aft of elastic axis, Xj/R 

offset of tension center aft of elastic axis, x^R 
(at the tip, XC should be set nearly equal Xl) 

2 2 

polar radius of gyration about elastic axis, kp/ R 

section moment of inertia about elastic axis, 

(slug-ft or kg-m) 

torsional st'.ffness, GJ (ib-ft^ or N-m^) 
blade twist relative .75R» (deg) 
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NLRTR 


OFSTLL 


OPYAW 

TAU(3) 

ADELAY 

AHAXNS 

psiasO) 

ALFDS(3) 

ALFRB(3) 

CLDSP 

CDDSP 

CMDSP 


Stall Bodel 

integer paraneter defining stall aodel 

0 no stall 

1 static stall 

2 McCroskey stall delay 

3 HcCroskey stall delay with dynamic stall 
vortex loads 

4 Boeing stall delay 

5 Boeing stall delay with dynamic stall 

vortex loads 

(the stall delay can be suppressed by setting TAU*0.) 

integer parameter defining jraved flow corrections 

0 both yawed flow and radial drag included 

1 no yawed flow (cosJV. = 1.) 

2 no radial drag (F “0.) 

3 neither yawed flow nor radial drag included 

stall delay time constants for lift, drag, and moment: 
(calculated if LT 0.) 

maximum angle of attack increment due to stall delay, 

•^max delay 

angle of attadc in linear range fcr no stall model, (deg) 

dynamic stall vortex load rise and fall time (azimuth increment) 
for lift, dra^, and moment: (deg) 

dynamic stall angle of attack for lift, drag, and moment: 

*^ds 

stall recovery angle of attadc for lift, drag, and moment: 
o<re (d««) 

maximum peak dynamic stall vortex induced lift coefficient: 
^‘^Sds 

maximum peak dynamic stall vortex induced drag coefficient: 

maximum peak dynamic stall vortex induced moment coefficient: 
^‘Hls 
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KHLNDA 

KFLMDA 

FXUtDA 

FnillM 

FMLHDA 

FACTWU 

KINTH 

KINTF 

KIMWB 

KINTHT 

KINTV'T 

HINGE 

NCXDLB 

NCCLT 

NDNROT 

EFN0I£ 


NLRT3 

factor for hover induced velocity (typically l.l) 

factor kj for forward flight induced velocity (typically 1.2) 

factor f for linear inflow vairiatlon in forward flight 
(typically 1.5) 

factor f for linear inflow variation in forward flight 
(typically 1.) 

factor f on linear inflow variation due to hub moment 
(typically 1.) 

factor introducing lag in C„, . and used to calculate 

induced velocity (typically .5) * ^ 


factor for hover interference velocity at other rotor 
(K. 2 J or H.^ 2 ) 


factor for forward flight interference velocity at other 
rotor (Kgi or K^ 2 ^ 

(lineeur variation between KINTH at m = O.O 5 and KINTF at 
^ * 0.10 is used) 


factor for rotor- induced interference velocity at 
wing-body, Ky 

factor for rotor- induced interference velocity at 
horizontal tail, 

factor for rotor- induced interference velocity at 
vertical tail, 

(Ky, Kjj, Ky equal fraction of fully-developed wake times 
maxinun fraction surface in wake) 


integer parameter specifying blade mode type 

0 hinged 

1 cantilever 

2 articulated (flap and lag nodes only) 

nunber of collocation functions for bending mode calculations 
(total flap and lag, alternatii^); maxinum 20 

number of collocation functions for torsion node calculations; 
maximum 10 

integer paranetert NE 0 to calculate nonrotating bending 
frequencies 

criterion on change of collective pitch to update blade 
modes, (deg) 
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NLRTR 


HASST 

XIT 

HBUEE 

EFLAF 

ELAG 

KFLAP 

KLAG 

RCPLS 

TSFRBG 

RCPL 

MDra 

WTIH 

FTO 

FTC 

FTR 

KTO 

KTC 

KTR 

KPIN 


PHI PH 
PHI PL 
RPB 
RPH 
XPH 


tip Bass (slug or kg)i the tip mass can also be included 
directly in the section nass distribution 

offset of tip Base center of gravity aft of elastic 
axis I x^/R 

blade Bass (slug or kg); if LB 0., integral of section nass 
used (with Bass included at r ■ 0. to account for the hub mass) 

flap hinge offset e^R (extent of rigid hub for cantilver blade) 

lag hinge offset e^/R (extent of rigid hub for cantilver blade) 

flap hinge spring (ft-lb/rad or B-N/rad) 

lag hinge spring (ft-lb/rad or n-N/rad) 

h.‘nge spring paraBster, 

S 

hinge spring paraBCter, 6^^ 

(hinge spring pitch angle is 

structural coupling paraBeterl^( effective pitch angle 1^0 
used to calculate blade bending nodes; nomally 1.) 

integer paraneter specifying twist inboard of r^^: EQ 1 for 

no pitch bearing 

integer paraneter defining control system stiffness input: 

1 for , 2 for 

control system frequency 0^(psr rev, at normal tip speed VTIPN) 
collective 

cyclic 

reactionless 

control system stiffness (ft-lb/rad or m-N/rad) 
colloctive 

cyclic 

reactionless 

integer parameter defining pitch/bending coupling input: 

1 for input, 2 for calculated (negative to suppress cosine 
factors in and Kp^) 

root geometry to calculate pitc)\/bending coupling (KPIN = 2 or - 2 ) 
pitch horn cant angle, (deg) 

pitch link cant angle, 

pitch bearing radial location, ipg/R 

pitch horn radial location, Tpj^R 

pitch horn length, *pjj/R 


-i4- 



atankp(hbm) 

DEL3G 

RFA 

ZFA 

XFA 

CONE 

DROOP 

SWEEP 

FEROOP 

FSWEEP 


f)RlC5INAL PAGE K ^ 

/ -1 , X QF POOR QUALITYi 

pltch/bendlng coupling tan Kp. (deg), for pitch 
horn level (KPIN » 1 or -l) ^ 

pitcV§i®^l coupling tan'^Kpp (deg), for pitch horn 
level 

feathering axis radial location, fp^/R 

giitibal undersling, Zp^/R 

torque offset, *p^/R 

precone angle positive up 

droop angle (^®s) ~ positive down from 

precone 

sweep angle ^fa^ (deg) at = 0 , positive aft 

feathering axis droop angle SFA 4 (*^®e)» positive down 
from precone 

feathering axis sweep angle ^fa^ (deg), positive aft 
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Namelist HLHAKE 


FACTUN 

OPVXVY 

KNW 

KRW 

KFW 

KIM 


RRU 

FRU 

PRU 

FNU 

DVS 

DLS 

CORB(5) 


0KX)RE(2) 


OPNWS(2) 


LHN 

OPHN 

OPRTS 


factor Introducing lag In bound circulation used to calculate 
Induced velocity 

Integer pauranetert EQ 0 to suppress x and y components of 
induced velocity calculated at the rotors 

extent of near wake, 

extent cf rolling up wake, 

extent of far wake and tip vortices, 

extent of far wake and tip vortices for points off rotor 
disk, Kjjj, 

(age 4 “ KAH' ; all K GE l) 

Initial ladlal station of wa’ce rollup, 

Initial tip vortex fraction of P ^ for rollup, fp„ 

max nu 

extent of rollup In wake age, (deg) 

tip vortex fraction of for near wake, f^jy 

sheet edge test parameter d^^{ LT 0. to suppress test 

lifting surface correction parameter d^^^; LT 0. to suppress 
correction 

vortex core radii r ^R 

( 1 ) tip vortices 

! 2j hurst tip vortices 
3; tip vortices In far wake off rotor 
(4 j trailed lines (LT 0. for default = s/2) 

(5) shed llr»s (LT 0, for default = t/2) 

Integer parameter specifying vortex core type: 0 for 

distributed vortlclty, 1 for concentrated vortlclty 
(l) tip vortices 

( 2 ) Inboard wake 

Integer parameter controlling action when Inflow and 

circulation points coincide In near wake ( 4 “ O) 

sheets are being used: 0 to use two sheets, 1 to use lines, 

2 to use single sheet 

(l) shed wake 
( 2 ) trailed wake 

number of spirals of far wake for axlsymmetrlc case, Lj^ 

Integer parameter: EQ 0 for axlsymmetrlc wake geometry 

Integer parameter: NE 0 to Include rotation matrices 

(R^, etc.) In Influence coefficients 


- 46 - 



.. * .iuii IS 

or FOUR Qb'AUTja 


NLtfAKE 

WKM0DL(13) Integer parameter defining wake model* 0 to omit 

element, 1 for line segment with stepped circulation 
distribution, 2 for line segment with linear circulation 
distribution, 3 for vortex sheet element 

(1) tip vortices (stepped line or linear line) 

(2) near wake shed vorticity 
p) near wake trailed vorticity 
p) rolling up wake shed vorticity 

(5) rolling up wake trailed vorticity 

(6) far wake shed vorticity 
(?) far wadte trailed vorticity 

(8) far waike (off rotor) shed vorticity 

(9) far wake (off rotor) trailed vorticity 

(10) bound vortices (no sheet model) 

(11) axlsymmetrical wake axial vorticity (no line model) 

(12) axlsymmetrical wake shed vorticity (no line model) 

(13) axlsymmetrical wake ring vorticity (no line model) 



NV.» gw 


MRG 

NG(MRG) 

MRL 

NL(MRL) 

0PWKBP(3) 


number of circulation points for near wake; LE MRA 

circulation points, identified by aerodynamic segment 

number: nQj^ for i = 1 to MRG (corresponding rj^ must be 

between r and 1.) 

root 

number of Inflow points; LE MRA 

points at which the induced velocity is calculated, 
identified by aerodynamic segment manberi n^. for 
1 = 1 to MRL 

integer parameter controlling blade position model for wake 
analysis 

(1) EQ 0 to suppress inplane motion 

(2) EQ 0 to suppress all harmonics except mean 

(3) EQ 0 for linear from r = ^root'^^ to r = 1. 
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NLUAKE 


VELB 

DPHIB 

IBV 

QDEBUG 


KRVG 

OHIWG 


FWGT(2) 

FWGSI(2) 

FHGSO(2) 

KWGT(4) 

KWGSI(4) 

KWGS0(4) 


core burst propagation rate, 

core burst age increment, (<ieg) 

core b\irst test parameter LT 0. to suppress bursting 

velocity criterion for debug print: print if 

\T.ic/r \ > (iEEBUG 


Prescribed wake geometry 

extent of prescribed wake geometry, K_„_ (age 4* “ KA.'V); 
maximum 144 

integer parameter defining prescribed wake geometry model 
1 from Kj=fj\ , , Input K^, input 

option #1, without Interference velocity in V 

from Input Kj^, K^, K^, 

Landgrebe prescribed wake geometry 

from 

from n___ 
max 

from V 

from V without interference 
Kocurek and Tangier prescribed wake geometry 
from C„ 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 


from rr 


max 


from \ 

from X without Interference 


Factors f. and f_ for prescribed wake geometry 
*tip vortex 

Inside sheet edge 

outside sheet edge 

Constants K. , K_, K^, for prescribed wake geometry 
tip vortex 

Inside sheet edge 

outside sheet edge 
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NLWAKE 


Free wake geometry 

KFWG extent of free wake geometry distortion calculation, K,™ 

(age suggest (.4 /m,)MPSI; maximum 96, 

multiple MPSI 

OPFWG Integer parameter defining free wake geoaetry model 

1 Scully free wake geometry 

2 option #1, without interfere nee velocity 

I'lERWG number of wake geometry iterations; suggest 2 or 3 

FACTMG factor introducing lag in distortion calculation to 

improve convergence; suggest O.5 

RTWG(2) radial station r/H of trailed vortlclty 

(l| inside sheet edge 

(2) outside sheet edge, or trailed line (suggest ,4) 

WGM0DL(2) integer parameter defining wake model: 0 to omit, 1 for 

line segment, 2 for sheet element 

(ll Iriboarxi trailed wake elements 
(2) shed wake elements 

CCREWG(4) vortex core radii r 

(1^ tip vortices 

(2) burst tip vortices (LE 0. for default = 
unburst value) 

(3) inboard trailed lines (LE 0. for default = 
i(RTWG(?) - RTMG(l)) ) 

(4) shed lines (LE 0. for default = 0.4A.V) 


MRVBWG 

LDMWG 

NDMWG(MPSl) 

DQWG(2) 


number of wake revolutions \:ised below point where induced 
velocity is being calculated; suggest 2 

integer parameter general update every increment 

in boundary age; suggest l80°/6Mi^ 

Integer parameter n^('V.)t boundary update every n^^ 

Increment in age, function of 3 = 1 to !^I; 

suggest fore and aft, and 45°/&H> on sides 

incremental velocity criteria; suggest 0.04 to O.OSX^ 

(1) near wake elements defined by 

> DQ¥G(l) 

(2) integrate bound vortex line in time over 
if |Aq| > DQWG(2) 


-^^original page is 

OF POOR QUALITY 



NLUAKE 


HVGDB(2) 


QWGIS 


Integer parameters controlling debug level 3 print 
of wake geometry distortion 

(1) IFRt print distortion before general 
update every IFR*a 1*} EQ 0 to suppress 

( 2 ) INPSi print distortion after each 
iteration every INFS*4&V> BQ 0 to 
suppress} last iteration printed in full 

parameter controlling debug level 3 printi induced velocity 
contribution of wake element printed if )A*q) > QVGIB; 
suggest 0.5X^ to 
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Namelist NLBODY 


TITLE(20) 

WEIGHT 

IXX 

lYY 

IZZ 

IXY 

IXZ 

lYZ 

TRATIO 

CONFIG 


ASHAFT(2) 

AGANT(2) 

ATILT 

HMAST 

DPSI21 

CANTHT 

CANTVT 


title for airframe and drive train data (80 characters) 

aircraft gross weight including rotcrs (ih or kg) 

aircraft moments of inertia including rotors (slug-ft'^ or kg-m"^) 
I 

XX 

I 

yy 

I 

zz 

I 

xy 

I 

xz 


yz 


ratio of rotor 
speed, 


M2 rotational speed t rotor #1 rotational 
(transmission gear ratio rj;^/rj2) 


integer pauameter specifying helicopter configuration 

0 for one rotor 

1 for single main rotor and tail rotor (rotor M2 
is the tall rotor) 

2 for tandem main rotors (rotor M2 is the rear rotor) 

3 for tilting pro pro tor aircraft (rotor M2 is the 
left rotor) 


shaft angle of attack (deg), positive rearward 

(1) rotor M\ 

(2) rotor M2 

shaft cant ang'' ■' (deg); positive to right for main 
rotor; posltiv<< upward for tail rotor; positive inward in 
helicopter mode for tilt rotor 

(1) rotor #1 

(2) rotor M2 


nacelle tilt angle o<.p (deg), for tilting proprotor 
configuration only; 0. for airplane mode, 90. for 
helicopter mode 

rotor mast length from pivot to hub (ft or m), for 
tilting proprotor configuration only 

(deg); rotor M2 azimuth amgle when rotor Ml 
angle * 0; must be 0, if ^ 1. 

horizontal tail cant angle 4^ (deg), positive to left 

vertical tail cant angle (deg), positive to right 


azimuth 
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NLBODY 


FSCG 

BLCG 

VLCG 

FSRl 

BLRl 

WLRl 

FSR2 

BLR 2 

WLR2 

FSMB 

BLtfB 

VIMB 

FSHT 

BLHT 

WLHT 

FSVT 

BLVT 

WLVT 

FSOFF 

BLOFF 

WLOFF 


location (fuselage station, butt line, and waterline) of 
aircraft components relative to a body fixed reference 
system having an arbitrary orientation and origin; fuselage 
station (FS) positive aft, butt line (BL) positive to right, 
and waterline (WL) positive up (ft or m) 

aircraft center of gravity location 


rotor #1 hub location (rl^t nacelle pivot location for 
tilting proprotor configuration) 


rotor #2 hub location 


wing-body center of action 


horizontal tall center of action 


vertical tall center of action 


^Int off rotor disk (for Induced velocity calculation) 


ClfrRLZ(ll) control Inputs (deg) for all stldcs centered (vp • 0)i 


-0 - ( 


0 nc Is 
rotor #1 




& 


Ic Is 
rotor #2 






aircraft 
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NLBCDY 


KOCFE 

KCCFE 

KSCFE 

KPCFE 

PCCFE 

PSCFE 

PPCFE 

KFOCFE 

KROCFE 

KFCCFE 

KRCCFB 

KFSCFE 

KRSCFE 

KFPCFE 

KRPCFB 

PFCCFE 

PRCCFE 

PFPCFE 

PRPCFE 

KFCFE 

KTCFE 

KACFE 

KECFE 

KRCP'E 


description of control system (for K parameters 

are gains (deg per stick deflection) parameters are 
swashplate azimuth lead angles (deg; 

one rotor, single main rotor and tail rotor, tilting proprotor 
configurations 

Kq, collective stick to collective pitch 

lateral cyclic stick to cyclic or differential 
collective pitch 

Kg, longitudinal cyclic stick to cyclic pitch 

K , pedal to tail rotor collective or dlffe-ential 
cyclic pitch 

, lateral cyclic stick to cyclic pitch (one rotor, 
or single main rotor and tail rotor configurations) 

, longitudinal cyclic stick to cyclic pitch 
s 

A.V , pedal to differential cyclic pitch (tilting 
propiotor configuration only) 

tandem main rotor configuration 

KpQ, collective stick to front collective pitch 

KpQ, collective stick to reair collective pitch 

Kjpq, lateral cyclic stick to front cyclic pitch 

KpQ, lateral cyclic stldt to rear cyclic p'tch 

Kpg, longitudinal cyclic stick to front collective pitch 

Kpg, longitudinal cyclic stick to rear collective pitch 

Kpp, pedal to front cyclic pitch 

Kpp, pedal to rear cyclic pitch 

lateral cyclic stick to front cyclic pitch 

^^RC’ cyclic stick to rear cyclic pitch 

AVpp» pedal to front cyclic pitch 

APpp, pedal to rear cyclic pitch 

aircraft controls (all configurations) 

Kp, collective stick to flaperon 

K^, collective stick to throttle 

K^, lateral cyclic stick to ailerons 

K , longitudinal cyclic stick to elevator 

Kj., pedal to rudder 
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NIBODY 


NBf 

QHASS(NE») 

QFR£(i(NBM) 

QDAHP(MEM} 

QDMIPA(HEM) 

QCmL(4,NEM) 

I10FSYK(MB«) 

ZETAR1(3*MEM) 

ZBTAR2(3,NEM) 

GAMAH1(3,1«M) 

GAMAR2(3*NEM) 


nvnber of alrfrane sodas for which data supplied t 
aaxiaus 10 

geoeralized aass including rotors (slug or kg) 
generalised frequency (Hz) 
structural daaplng g^ 

^e^ynanjc dasplng « i(Q,/i 5 V^)/^(q 5 j^V) 

control derivative F- # “o i(Qw/'l\V^)/i^ for 

Se» Sa* 6r tk/n^ 

integer vector designating type of aodet GT 0 for 
synnetric, LT 0 for antisyaaetric; only required for 
flutter analysis wi^ OFSYMN MB 0 

linear node sh pe at rotor #1 huh (ft/ft or n/n) 

linear node shape at rotor #2 hub (ft/ft or n/n) 

angular node shape at rotor #1 hub (rad/ft or rad/«) 
angular node shape iTj^ at rotor 02 hub (rad/ft or rad/m) 



pitcb/nast-bending coupling (rad/ft or rad/m) 

kfhci(neh) 

KMCk “ - 


for rotor 01 

KFNSl(HEM) 

KMSk “ - 

^®ia/^qsk 

for rotor #1 

KFHC2(NEM) 

KMCk - - 


for rotor #2 

KFHS2(MBN) 

KMSk “ • 


for rotor #2 
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Aircraft aerodynaaic cteracterlstics 
Wing-body 


LFTAW 

L^/q 

(ft^rad or 

LFTFW 

Lif/q 

(ft^/rad or mV^ad) 

LFTDW 


(ft^rad or m^/rad) 

AMAXW 

*^max 

(deg) 

IWB 


(deg) 

meow 

"WB “ 

(ft^ or m^) 

IRGVW 

^vert„ - 

(ft^ or m^) 

IRGIW 

n«S^ - { )■' 

(ft^ or m^) 

LRGFW 

°0h/i 

(ft^rad or m^/rad) 

IRGDW 


(ft^/rad or m^/rad) 

MOMOW 

M(/q 

(ft^ or m^) 

MOMAV 

M^/q 

(ft^/rad or m^/rad) 

MOMFW 

Mjf/q 

(ft^/rad or mV^ad) 

MOMDH 

w%A 

(ft^/rad or m^/rad) 

SIIEB 


(ft^/rad or m^rad) 

SIDEP 

VYp/q 

(ftV^ad or m^/rad' 

SIIER 

VY^q 

(ft^/rad or m^/rad) 

ROLLS 

Nx^q 

(ft^/rad or ra^/rad) 

ROLLP 

vNxp/q 

(ft^/rad or m^/rad) 

ROLLR 

VNx/q 

(ftV^ad or mVrad) 

ROLLA 

Nx^/q 

(ft^/rad or m^/rad) 

YAWS 

YAWP 

vNzp/q 

(ft^/rad or m^rad) 
(ftVrad or mV^ad) 

YAWR 

vNz/q 

(ftVrad or m^rad) 

YAWA 


(ft^/rad or m^/rad) 


Horizontal tall 


LFTAH 

Lo</q 

(ft^/rad or mV^^d) 

LFTEH 


(ft^rad or m^rad) 

AMAXH 

^max 

(deg) 

IHT 


(deg) 
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LFTAV 

LFTRV 

ANAXV 

IVT 

pmit 

LHTAIL 

HVTAIL 

OFTIirr 


Vertical tall 



(ft^/xBd or m^/r»d) 
(ft^rad or a^rad) 
(deg) 

(dag) 


Alrfraae interference 

horisontal tail length 5^^ for € (ft or ■) 

vertical tail height for positive up (ft or a) 

integer parameter controlling alrfxane/ tail aerodynamic 
interference! BQ 0 to suppreae (c “ 0 and ^-0) 
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EMGFOS 


lENG 

KMASTl 

KMAST2 

KICS 

KENG 

GSE 

GSI 

KEDAMP 

THRTLC 


KPGOVB 

KPGOVl 

KPG0V2 

KIGOVE 

KIGOVl 

KIGCV2 


Engine and drlTe train paxaaeters 

Integer paraaeter specifying drive train configiirationi 

0 one rotor 

1 asynetric, engine Isy rotor #1 

2 asynnetric, engine 'by rotor ^ 

3 syuetric 

2 

engine rotational inertia ^for both 2 engines if 

symetric configuration (slug-^t*" or kg-n^) 

drive train spring constants (ft-lb/rad or a-N/rad) 
rotor #1 shaft, Kmi or 

rotor #2 shaft. Kn 

^2 2 
interconnect shaft, or r^ 

engine shaft, r^ 

engine shaft structural damping g ( degree of freedom) 

S 6 

interconnect shaft structural damping ( *Pj degree of 
freedom) 

engine damping factor K: typically 1.0 for turhoshaft 
engines, or 10. for Induction electric motors 

(dimensional), for both engines if syaaetric 
configuration; if the throttle variable is only used 
for the governor, just the products 

Kl 

must be correct (P = 

governor proportional feedback gains (sec) 
to throttle, Kp * - 

to rotor #1 collective, Kp = 

to rotor #2 collective, Kp = 

governor integral feedback gains 
to throttle, K^ = - 

to rotor #1 collective , = A0/ 

to rotor #2 collective, K_ = ^9/dV., 

X s 
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TIGOVB 

TlQOVl 

T1G0V2 

T200VB 

T200V1 

T2Q0V2 


Soveroor ti«e las (sec) 

throttle ^ 

rotor #1 

rotor #2 

governor tiae lag (®*c ) 

throttle 

rotor #1 

rotor #2 
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Naaellst HLLOAD 


MVIB 

fsvib(mvib) 

BLVIB(MVIB) 

wlvib(mvib) 

ZETAV(3,NBM,MVIB) 

HALOAD 

MHLOAD 

MRLOAD 

RLOAD( MHLOAD) 
MHARML 

NPOLAR 

NWKGMP(4) 

HWKGNF 

J'-fKGMP(MWKGMP) 


Alrfraae vibration 

number of stations for airframe vibration 
calculation and print; maximum 10; LB 0 to 
suppress 

airframe location for vibration calculation (ft or m) 
fusel£ige station 

butt line 

waterline 

linear node shape airframe vibration 

stations (ft/ft or m/m) 


Integer parameter controlling print of motion and 
aerodyneunlcst BQ 0 to suppress; LT 0 for only plots 

Integer parameter controlling print of hub and 
control loads: BCl 0 to suppress 

number of radial stations for blade section load 
calculation and print; maximum 20; LB 0 to suppress 

blade radial stations r/R for section loads 

ntimber of harmonics In loads analysis; maximum 30 ; 

LT 0 for no harmonic analysis; suggest about MPSl/3 

Integer parameter n for polar plots: symbol printed 

every n-th step 

Integer parameter controlling wadte geometry printer 
plot; BQ 0 to suppress 

( 1 ) top view 

( 2) side view 
back view 

(4) axial convection 

number of azimuth stations at which wake geometry 
plotted; maximum 8; LE 0 for no plots 

azimuth stations at which wake geometry plotted 
(S>- jAS>) 
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NPLOT(75) 


Integer parameter controlling printer-plots of motion 
and aerodynamics I 0 for no plot, 1 for time history 
plot, 2 for polar plot, 3 for both (only time history 
available for 1-4 and 68-75) 


( 1 ) bending motion 

(2) torsion motion 

(3) maximum circulation 

(4) V off rotor 

(5) << 

(6) M 

(7) A 

(8 c, 

(9) Cd 


Ul] 

:i2 


^dradlal 

P 


(13) up 

(14) UT 

(15; UR 

(1 6 ) u 

(17) © 

( 18 ) 

(19) lag 

(20) flap 

\lh 

(23 


^eff 


drag 

moment 


24) lift 

drag 


moment 


interference \ 


(25 
26 
(27 
(28 
(29 

30 

31 

32 

(33 ) Uq 

(35) wg 

(36) L/c 

(37) D/c 

(38) m/c 

(39) Dr/c 

(40) Fx/c 

(41) Fj,/c 

(42) Fa/c - CJ^ 

(43 Ma/c 
(44) 
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( 68 ) 

69 ) 

(70) 

(71) 

(72) 

(73) 

(74) 

(75) 


not used 
not used 
not used 



» 

« 

« 

« 

* 

* 

« 

« 

* 


not used 
not used 
not used 
P * 



rotating frame root loads 
nonrotating frame hub loads 
rotating frame root loads 
nonrotating frame hub loads 
section loads, shaft ajces 
section loads, principal axes 
section loads, shaft axes 
section loads, principal axes 


« 

* 


* 


^dimensional quantities 

for polar plots, last digit of Integer part of 
(value/increment) is printed, if it is a multiple 
of NFOLAR; the plot increment is defined as follows 
.01 plots 27-35 

.1 plots 6, 8-16, 24-26, 36-51 

1. plots 5. 7, 17-23, 52-61 

10. plots 62-67 
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KFATIG 

SBNDUR(18) 

CMAT(18) 

EXMAT(18) 


parameter K in fatigue damage calculation; suggest 
3 or 4 

endurance limit Sg (dimensional force or moment) 
material constant C 


material exponent M 


rotating frame root loads 
inplane shear f 
axial shear f * 
vertical shear f^ 
flap moment m 
lag moment m 
control moment m 
nonrotating frame hub loads 
drag force H 
side force Y 
thrust T 
roll moment M 
pitch moment M 
torque Q ^ 

section loads (principal axes) 
chord shear f^ 
axial shear f 
normal shear f 
flatwise moment m 

X 

edgewise moment m^ 
torsion moment m. 


the S-N curve is approximated by N = C / (S/Sg- l)^ 

use S„ LT 0. or G LT 0. to suppress damage fraction 
calculation; use M EQ 0. to suppress equivalent 
peak-to-peak load calculation as well 
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Far field rotational noise 


MNDISE 

RANGE(MNOISE) 

ELVATN(MMOISE) 

AZMUTH(MHOISE) 


num^ber of microphones; maximum 10; LB 0 for no 
noise analysis 

microphone range relative hub (ft or m) 

microphone elevation relative hub (deg), positive 
above rotor disk 

microphone azimuth relative hub (deg) , defined as 
for rotor azimuth 


MHARMN(3) 


MTIMEK(3) 


AXS(MRA) 

OPNOIS(4) 


number of harmonics 

(l) in noise calculation; maximum 500 
(2; in aerodyncunic load harmonic analysis 
(suggest MPSl/3) 

(3) in print of noise (LE 0 for no print) 

number of time steps (LE 0 to suppress) 

(1) in period of noise calculation; maximum 500 

(2) increment in noise print 

(3) increment in noise plot 

blade cross section area Aj-g/c*^ at aerodynamic 
segments, for thickness noise calculation( typically 
0.685 times thidcness ratio) 

integer parameter controlling noise calculation: 

0 to suppress, 1 for impulsive chordwise loading, 

2 for distributed chordwise loading 
(l) lift noise 
(2 j drag noise 
(3) radial force noise 
(4; thidcness noise 
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Namelist NLFLUT 


OPFLOW 

OPSYMM 

OFF DAN 
NBLEFL 
MPSIFC 


NINTFG 

MFSICC 

DALPKA 

rWACH 

OFUSLD 

DELTA 

OFRINT 

OFGRND 

KASGB 

OFSAS 

KCSAS 

KSSAS 

TGSAS 

TSSAS 


integer parameter specifying analysis type: LT 0 for 

constant coefficient approximation; BQ 0 for axial flow; 
GT 0 for periodic coefficients 


integer parauneter: NE 0 for symmetric and antisymmetric 

analyses (tilting proprotor configviration only) 

integer parameter; EQ 0 to suppress flight dynamics analysis 

integer parameter: EQ 1 for independent rotor blade analysis 


number of azimuth steps in period for nonaxial flow, periodic 
coefficient analysis (OFFLOW GT O); = 360/(N,^,M) for 

odd number of blades, = 720/(N . ,m) for even number of 

blades 


integer parameter specifying numerical integration option 
for periodic coefficient analysis (OFFLOW GT O) : 1 for 

modified trapezoidal method, 2 for Runge-Kutta method 

number of azimuth stations (per revolution) in evaluation 
of average coefficients for constant coefficient approximation 
(OFFLOW LT 0); 360°/M 

angle of attack increment .Aot (deg) for calculation of 
Cj^, c^, and c^ derivatives in aerodynamic coefficients 

Mach number increment AJ1 /m for calculation of Cj^, c^, and 
c^ derivai.ives in aerodynamic coefficients 

integer parameter controlling use of unsteady lift and moment 
in flutter analysis: 0 to suppress; 1 to include; 2 for 

zero in stall (15°'< < 165°) 

control and motion increment for aircraft stability derivative 
calculation (dimensionless) 

integer parameter: EQ 0 to suppress rotor/body aerodynamic 

interference in flutter analysis 

integer parameter controlling ground effect analysts; EQ 0 
for out of ground effect, NE 0 for in ground effect 

factor for antisymmetric ground effect; 0. to suppress, 1.0 
for unstable roll moment due to ground effect (tilting proprotor 
configuration only) 

Integer parameter controlling use of SAS: EQ 0 to suppress 

lateral SAS gain = - (deg/deg) 

longitudinal SAS gain (deg/deg) 

lateral SAS lead time (sec) 
longitudinal SAS lead time X„ (sec) 
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NLFLUT 


0PTORS(2) 

D0F(80) 


C0H(26) 


GUS(3) 


Integer parameter i EQ 0 for rigid pitch model (infinite 
control system stiffness, no degree of freedom) 

(1) rotor #1 

V2) rotor #2 


integer vector defining degrees of freedom for flutter 
analysis; 0 If not used, 1 If used, 2 If quaslstatlc variable; 
order I 


rotor #1 
rotor #2 


At) tn cn 


iVSCti CV) 

^ d|C ®ii 


**»/» 


ftC) 

bending 

(15) 






•C. 


pitch/ tors ion 

(9) 


f&fc Vc 

gimbal rotor Inflow 
teeter si>eed 


alri.ame 

rigid body flexible engine governor 
body ( 10) speed 


integer vector defining control variables, 0 if not used; 
order* 

rotor #1 ©o Sic ©•% •— 

rotor #2 ©« Sic ••• 

pitch (8) 

airframe Sc 

pilot S* i j 

Integer vector defining gust components, 0 If not used; 
order* u^, Vq, w^ 

for a two-bladed rotor, Is replaced by 

tnere are rotor pitch control variables; except for a 

two-bladed rotor, which has the 4 variables 9^^, 9^^^, 
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NLFLUT 


ANTYPE(4) 


NSYSAN 

NSTEP 

NFREQ 

FREQ(NFRBQ) 


Integer parameter specifying tasks in linear system 
analysis, 0 to suppress 

(1) eigenanalysis 

(2) transfer function printer-plot 
(3; time iilstory printer- plot 

(4) rms gust response 

Eigen 2 malysis 

calculation control: 0 for eigenvalues, 1 for eigenvalues 

and eigenvectors; 10 or 11 for zeros as well 

static response calculated if NE 0 

number of frequencies for which frequency response 
calculated; LE 0 to suppress; maximum 100 

vector of frequencies (per rev) 


NBPLOT 

NXPLT 

WPLT 

NAMEXP( NXPLT) 

namevp(nvplt) 

NDPLT 

NFOPLT 

NFIPLT 

MSPLT 


Transfer function printer-plot 

calculation met*«)d: if 1, from matrices; if 2, from 

poles amd zeros 

maaber of degrees of freedom to be plotted; maximum 80 

number of controls to be plotted; maximum 29 

vector of variable names to be plotted (inconsistent 
names ignored ) 

vector of contr'l names to be plotted (inconsistent 
names ignored) 

frequency steps per decade 

exponent (base lO) of beginning frequency 

exponent (base lO) of end frequency 

(maximum NF= (NFIPLT - NFOPLT) * NDPLT + 1 = 150 

maignitude plot scale: if 1, plot relative maximum value; 

if 2, plot relative 10.**K; if 3» plot relative 10. 
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NLFLUT 


Time history prlnter-ilot 


NTPLOT 

control input type: 1 for step, 2 for impulse, 3 for 

cosine impulse, 4 for sine doublet, 5 for square impulse, 
6 for square doublet 

PERPLT 

period T for impulse or doublet (sec) 

IFPLT 

time step (sec) 

TMXPLT 

maximum time (sec); maximum NXPLT»N'/PLT*TMXPLT/DTPLT= 7200 

NXPLT 

number of degrees of freedom to be plotted; maximum 80 

IT/PLT 

number of controls to be plotted; maximum 29 

NAMEXP( NXPLT) 

vector of variable names to be plotted (inconsistent 
names Ignored) 

namevp(nvplt) 

vector of control names to be plotted (inconsistent 
names ignored) 


Rms gust response 

lgust(mg) 

real vector of gust correlation lengths: GT 0., dimensional 

length L ('Cg = L/2V); EQ 0., set L = 400.; LT 0., 
magnitude is^dlmensionless correlation time (frequency 

l- - 

mgust(mg) 

real vector of gust component relative magnitudes 


MG = number of gust components; maximum 3 

namexa(macc) 

vector of names of degrees of freedom for which 
acceleration calculated; last 3 eqioal AGGB for body 
axis acceleration (all 3 or none) (inconsistent names 
ignored) 

FREQA(MACC) 

vector of acceleration breadc frequencies (Hz); 2/rev 
use<’ if LT 0.; in seime order as NAMSXA 

MAGC 

number of accelerometers; LE 0 for none; maximtun 83 

FSACC 

location of point at which body axis acceleration 
calculated (ft or m) 

fusel 3ge station 

BLACG 

butt line 

WLACC 

waterline 

zbtacc(3,nem) 

linear mode shape at point where body axis acceleration 

calculated 

NAMEXR(3) 

names of state vector; assumed 

that Sis* ®ls Immediately ( Inconsistent 

names ignored; 
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NLFL13T 


Vairlable names for llneaur system analysis 


Degrees of freedom 



IBl 

• 

• 

1B15 


?*»^x 

bending 

ITl 

e 



pitch/ torsion 

• 

1T9 




IBGC 

p<ic 


gimbal/ teeter 

IBGS 




PSIS 



rotor si>eed 

ILU 

Vik 


inflow 

lU 




ILY 




ZBl 

• 

• 

-A 

bending 

ZB15 


li'i 

•' 9n/i 


2T1 

e 

sli'eli'e " . 

pitch/ torsion 

• 




ZBGC 



gimbai/ teeter 

ZBGS 

(3e>& 



PSII 



rotor speed 

2LU 

Vu 


inflow 

2LX 




2LY 




PHIF 



rigid body 

THTF 

ep 



PSIF 

Vp 



XF 




YF 




ZF 





rotor #1 


1 

rotor #2 


airframe 


V 
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NLFLUT 


QFl 

• 


flexible body 

• 

• 

QFIO 




PSIE 


engine speed 

TGOV 


governor 

IGOV 




2G0V 




Control 

variables 



ICC 


rob 

3r #1 

ICIC 




ICIS 



i 

1C4 

• 




• 

1C& 



1 

t 

2C0 

®IS ^I^/X 

rob 

sr #2 

2C1C 




2C1S 




2C4 




CO 




EELF 


airframe 

i£LB 




EELA 


1 


QBLR 

ir 

1 


CT 

©^ 

1 

1 

lELO 

i. 

pilot 

DELC 

Sc 



EELS 




EELP 

Se 



EELT 

St 


r 


airframe 


V 
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NLFLUT 


Gust components 
UG 

VG V, 

KG 

the leading character (l or 2) Is 
depending on the helicopter configuration 
blank (left jiistifled) 

M or T 
F or R 

R or L (OPSYMM = O) 

S or A (OPSYMM f 0) 

For a two bladed rotor, BGC is replaced by BT 

For first order degrees of freedom, the only state is the 
velocity, hence It is the velocity that will be plotted 


For the rotor names, 
replaced as follows, 
CONFIG = 0 
1 
2 
3 
3 
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Naaelist NLSTAB 


NPRNTP integer parameter controlling performance print during 

stability derivative calculation: LE 0 to suppress 

NPRNTL integer parameter controlling loads print during stability 

derivative calculation: LE 0 to suppress 

ITERS number of '^adce influence coefficient/motion and forces 

iterations 

OPLMDA integer parameter controlling iiKiuced velocity calculation: 

if 0, update Influence coefficients and inflow; if 1, suppress 
influence coefficient update; if 2, si ppress inflow update 
(and influence coefficient update) 

lELTA control and motion increment for stability derivative 

calculation (dimensionless) 

D0F(?) integer vector defining degrees of freedom, 0 if not used; 

order: 4^^, dp. Xp, yp, Zp, 

C0N(16) Integer vector defining control variables, 0 if not used; 
order: 

rotor #1 

rotor #2 3o Oic 
airframe ^ ®t. 

pilot 


GUS(3) integer vector defining gust components, 0 if not used; 

order: Uj., v^, w^ 

CPPRNT(4) integer parameters controlling stability derivative print, 
EQ 0 to suppress: 

(1) rotor coefficient form, dimensionless 

(2) rotor coefficient form, dimensional 

(3) stability derivative form, dimensi ""'ess 

(4) stability derivative form, dlmensii 

KCSAS lateral SAS gain, (deg/deg) 

KSSAS longitudinal SAS gain, K = (deg/deg) 

TCSAS lateral SAS lead time (sec) 

TSSAS longitudinal SAS lead time (sec) 
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EQTYPE(12) 


ANTYFE(5) 


NSYSAN 

NSTEP 

NFREQ 

FREQ(NFREQ) 


.rO^STAB 


Integer parameter specifying equations to be 
analyzed, EQ 0 to suppress 

with I with SAS 

(1) complete 

(2) symmetric 

( 3 ) antisymmetric 
with ^ , without SAS 

(4) complete 

( 5 ) symmetric 

(6) antisymmetric 
without V , with SAS 

( 7 ^ complete 

(R> symmetric 

( 9 ) antisymmetric 

without without SAS 

( 10 ) complete 

( 11 ) symmetric 

( 12 ) antisymmetric 


integer parameter specifying tasks in linear system 
analysis, EQ 0 to suppress 

( 1 ) eigenanalysis 

( 2 ) traiisfer function printer-plot 

(3) time history printer-plot 

(4) rms gust response 

( 5 ) numerical integration of transient 

Eigenanalysis 

calculation control: 0 for eigenvalues, 1 for eigenvalues 

and eigenvectors; 10 or 11 for zeros as well 

static response calculated if NE 0 

number of frequencies for which frequency response 
calculated; LE 0 to suppress; maximum 100 

vector of frequencies (per rev) 
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NLSTAB 


Transfer function printer-plot 

NBPLOT calculation aethod: if 1, from matrices; if 2, 

from poles and zeros 

NAMEXP(NXPLT) vector of variable munes to be plotted (inconsistent 
names ignored) 

NAMEVP( NVFLT) vector of control names to be plotted (inconsistent 
names ignored) 

NXPLT number of degrees of freedom to be plotted; maximum 7 

NVPLT number of controls to be plotted; maximum 19 

NDPLT frequency steps per decade 

NFOPLT exponent (base lO) of beginning frequency 

NFIPLT exponent (base lO) of end frequency 

(maximum NF= (NFIPLT NFOPLT) * NDPLT + 1 = 15l) 

MSPLT magnitude plot scale: if 1, plot relative maximum 

value; if 2, plot relative 10**K; if 3t plot relative 10. 

Time history printer-plot 

NTPLOT control input type: 1 for step, 2 for impulse, 3 for 

cosine impulse, 4 for sine doublet, 5 for square impulse, 

6 for square doublet 

PERPLT period T for Impulse or doublet (sec) 

DTPLT time step (sec) 

TMXPLT maximum time (sec); maximum NXPLT*NVPLT*'IMXPLT/DTPLT = 7200 

NXPLT number of degrees of freedom to be plotted; maximum 7 

NVPLT number of controls to be plotted; maximum 19 

NAMEXP( NXPLT) vector of variable names to be plotted (inconsistent 

names ignored) 

NAMEVP( NVPLT) vector of control names to be plotted (inconsistent 

names ignored) 

Rms gust response 

LGUST(MG) real vector of gust correlition lengths: GT 0., 

dimensional length L (TTq == b/2V); EQ 0., set L = 400.; 

LT 0., magnitude is dimensionless correlation time 
(frequency U 

MGUST(MG) real vector of gust component relative magnitudes 

MG “ number of gust components, maximum 3 
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NLSTAB 


NAMSXA(MACC) 

vector nf names of degrees of freedom for which 
acceleration calculated; last 3 equal ACCB for body 
axis acceleration (all 3 or none) (inconsistent names 
Ignored) 

freqa(macc) 

vector of acceleration break frequencies (Hz); 2/rev 
used if LT 0.; same order as NAMEXA 

MACC 

number of accelerometers; LB 0 for none; maximum 10 

location of point at which body axis acceleration 
calculated (ft or m) 

FSACC 

fuselaige station 

BLACC 

butt line 

WLACC 

waterline 

Numerical integration of transient 

TSIEP 

time step in numerical integration (sec) 

TMAX 

maximum time in numerical integration (sec) 

NPRNTT 

integer parauneter nt transient print every n-th 
integration step; LE 0 to suppress 

OPPLOT 

integer parameter controlling printer plot of body 
motion: EQ 0 to suppress 

D0FPLT(21) 

integer vector designating variables to be plotted, 

EQ 0 if not plotted; order; 

^ ipy 4f Jp ^ % 'Uy % 

OPTHAN 

see namelist NLTRAN 


CTIME 

CMAG(5) 

GTIME 

GMAG(3) 

GDIST(2) 

VELG 

PSIG 

0PGUST(3) 
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NLSTAB 


Variable names 

for linear 

system analysis 

Degrees of freedom 


PHIF 


rigid body 

THTb 



PSIF 



XF 



YF 



ZF 



PSIS 


rotor speed 

Control 

variables 


ICO 


rotor #1 

IGIC 



IGIS 



2G0 


rotor ;f2 

2G1G 



2G1S 



lEL? 


aircraft 

DELE 



DELA 

u 


DELR 



GT 

9x 


mo 


pilot 

DELG 



DELS 



DELP 



BELT 



Gust components 


UG 



VG 



WG 
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NLSTAB 


For the rotor control names, the leading character (l or 2) 
is replaced as fellows, depending on the helicopter configuration 


CONFIG = 0 
1 
2 
3 


blank (left justified) 
M or T 
F or R 
R or L 


For first order degrees of freedom the only state is the 
velocity; hence it is the velocity that will be plotted 
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Namelist HLTOAN 


nRtGINAL’ PAGE Rf 

• ' M? on-iUTv 


NPRiN'TT 

NPRNTP 

NPRNTL 

NHSTRT 

TSTEP 

TMAX 

ITERT 

OPLMDA 


D0F(7) 


OFSAS 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

OPPLOT 

D0FPLT(2l) 


integer parameter n; transient/performance/loads print 
every n-th integration step; LE 0 to suppress 

integer parameter controlling j)erformance print: LE 0 to 
suppress 

integer parameter controlling loads print: LE 0 to suppress 

integer parameter n: restart file written only every n-th 

integration step; LE 0 to suppress 

time step in numerical integration (sec) 

maximum time in numerical integration (sec) 

number of wake influence coevficients/motlon and forces 
Iterations 

integer parameter controlling Induced velocity calculation: 
if 0, update influence coefficients and inflow; if 1, suppress 
influence coefficient update; if 2, suppress inflow update 
(and influence coefficient update) 

integer vector defining degrees of freedom in numerical 
integration; EQ 0 to suppress acceleration; 
order: vVp, x^., y^, Zp, 

integer parameter controlling use of SAS: EQ 0 to suppress 

lateral SAS gain, = - (deg/deg) 

longitudinal SAS gain, (deg/deg) 

lateral SAS lead time (sec) 

longitudinal SAS lead time f_ (sec) 

s 

integer parameter controlling printer plot of body motion: 

EQ 0 to suppress 

integer vector designating variables to be plotted; EQ 0 for 
not plotted; order: 
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NLTRAN 


OPTRAN 

CTIME 

CMAG(5y 


GTIME 

GMAG(3) 


GDIST(2) 

VELG 

PSIG 

OPGUST(3) 


Transient gust and control 

Integer parameter specifying transient option; 1 for 
control; 2 for uniform gust: 3 for convected gust 

period T for control (sec) 

control magnitude vp^ = ( Sq S^) (deg) 

defines cosine control transient with period T 
and magnitude vp^ 

period T for uniform gust (sec) 

gust maignitude g^ = (u^ w^) (ft/sec or m/sec) 

defines cosine uniform gust transient with 
period T and magnitude 

lengths for convected gvist (ft or m) 

(1) wavelength L 

(2) starting position 

gust convection velocity V (ft/sec or m/sec) 

D 

azimuth angle of convected gust wave front V (deg) 

s 

integer parameters defining convected gust model 

(1) EQ 0 to not use 

(2) rotor #1: 0 for gust at hub, 1 for over disk 

(3) rotor #2; 0 for gust at hub, 1 for over disk 

defines cosine convected gust transient with 

wavelength L and magnitude g ; for L = R the 

wave starts at edge of rotor disk, for L = 0. 

the wave starts at hub — assuming the aircraft 

center of gravity is directly below the hub; 

convected at rate V relative to moving aircraft 

if V is not used. It rate V relative to fixed 

frame if V is used ® 

a 
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ORIGINAL PAGE IS 
OF POOR OP \[.n V 

NLTRAN 


Transient gust and control subroutines 

The subroutine CONTRL calculates the transient control time 
history, C(t). The subroutine GUSTU calculates the uniform 
gust time history, G(t). The subroutine GUSTO calculates 
the convected gust wave shape, G(x ). The subroutines 
presently calculate a cosine-impulle gusts 

CONTRL C(t) = i(l - cos2TTt/T ) 

GUSTU G(t) =1(1 - cos2'irt/T) 

GUS"G G(x ) = id - cos2n(x -L )/L ) 

S so 

Other transients may be use^ by re^^acing these subroutines 
as required. 
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Namelist Inputs for Old Job (Restart) 


Namelist NLTRIM 

ANTYPE(3) 

OPREAD(IO) 

DEBUG(25) 

NHINTI 

Namelist NLFLUT 

ANTYFii(4) 

NSYSAN 

NAMEXR(3) 

Namelist NLSTAB 

0PPRNT(4; 

KCSAS 

KSSAS 

TCSAS 

TSSAS 

EQTYPEC12) 

ANTYPE(5) 

NSYSAN 

'"PGUST(3) 

Namelist NLTRAN 

NPRNTT 

NPRNTP 

NPRNTL 

NRSTRT 

TMAX 
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7 . NOTES ON PRINTED OUTPUT 


ORfGINAU PAGE IS 
OF POOH QliAfJTY 


This section presents notes on the printed output of the program i 
particularly regarding the units of the variables appearing in the output. 

Print of Performance (Program P ERF) 

Operating conditlonj 

a) motion; 1st number imensionless, 2nd number dimensional 

1) velocity - ft/sec or m/sec 

2 2 

2) dynamic pressure, q = Ib/ft or N/m 

3) weight, = lb or N 

4) body motion « deg/sec, ft/s'^c or m/sec 

5) *z = tt/sec or m/sec^ 

6) Vg = rpm 

b) body orientation and controls in deg 
Circulation convergence s 

a) tolerance, CG/S in form 

b) G/E = ratio error to tolerance ( ^ 1. if converged) 

Motion convergence: 

a) tolerance, BETA (etc) in deg 

b) BETA/B (etc) = ratio error to tolerance (^ 1. if converged) 
Airframe performance; section ^.2.6 

a) aerodynamic loads: dimensional 

b) components 

1) angles in deg 

2) loads, q dimensional 

3) induced velocity, total velocity dimensionless 
Gust velocity; dimensionless 

System power; 

a) dimensional (HP)i ntwiber in parentheses is percent total power 

b) climb power “ V^W 
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SysttA efficiency per'aneters: 

a) gross weight, VI “ lb or H 

b) drag-rotor * D/q-rotor = D 

L/D-rotor “ VI/D ^ 

c) drag-total - = ”tolalA' “A-*”**! ' “tota/H' ^ 

l/D-total - 

d) flgu::, of aarlt - « - 1 - 


Print of Rotor Loads (Program LCAlBl) 

Print aerodynamics v function r and H^) 

a) dimensionless quantities generally, angles in degrees 

b) induced velocity in nonrotatlrg shaft axes (\„, "“X-) 

'X y 'z 

c) interference induced velocity is that due to other rotor 

d) gust components i. ’city axes 
Force/c (dimensionless) > 


L/C - 5u = L/o^ 

I*/C - -V=, 

M/C - - M/o. 


SB/C - iU'^(c/c.^„)c4 


mean ^ ^radial 


D /c 
radiar mean 


F2/C = CT/S - F /c^ » d(C_yv}/dr 

% mean i 

FX/G = F /c 
' r mean 


Ma/C * M,/c 
' a' mean 


-VVan 

FBT/C - ?^/c,^„ 


-S/ - 



ORIGINAL PAGE IS 
OF POOR quality 


Forces (dlnenslonal) 


L ® section lift 

(Ib/ft or 

N/m) 

D = section drag 

(Ib/ft or 

N/m) 

M = section pitch moment 

(ft-lb/ft 

or m-N'/n) 

DR = section radial drag 

( Ib/ft or 

N/m) 

FZ = F = dT/dr 
z ' 

( Ib/ft or 

N/m) 

X 

11 

X 

(lb/ 'ft or 

N/m) 

MA = M 

a 

(ft-lb/ft 

or m-N/m) 

FR = F 

r 

(Ib/ft or 

N/m) 

FRT = F 

r 

( Ib/ft or 

N/m) 

Blade section power: section 5 •2.1 



CP/S = d(Cp/v-)/dr 




P = section power (HP/ft or HP/m) 


Print During Stability Derivative Calculation (Program STASH) 

a) Increment; 1st number diiaensionless, 2nd number dimensional 

b) motion and controls : 1st number dimensionless, 2nd number 

di' snsional 

1) angulsir velocity = deg/sec 

2) linear velociuy, gust velocity = ft/sec or m/sec 

3) = rpm 

4) Zp = ft/sec^ or m/sec^ 

5) controls = deg 

o; generalized forces: moments and forces in ^2C/v~a. form 

(rotor M parameters, body axes); torque in -<SG_/«-a 
form (rotor #1 peu^meters) 

Print of Stability Derivatives (Program STABP) 

Cp+ions ; 

a) rotor coefficient form, M X = 32C/v-a 

b) stability derivative form, X (acceleration) 

c) dimensionless or dimensional 
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b) subscripts 

acceleration (z) » rfR (FA) 

angular velocity = STL 

linear velocity = (FV) 

controls = 57 .3 

gust velocity = riR (FV) 


Print During Flight Dynaaics Numerical Integration (Program STABP 

a) controls in deg 

b) gust velocity: 1st number dimensionless, 2nd number 

dimensional 

c) aircraft motion: 1st number dimensionless, 2nd number 

dimensional 

1) displacement = deg, ft or m 

2) velocity = deg/sec, ft/sec or m/sec 

3) acceleration = deg/sec , g 

4) inertial axes = deg/sec, g 


Print Transient Solution (Prograun IBANPy 

a) controls in leg 

b) gust velocity dimensional 

c) aircraft motion: 1st number dimensionless, 2nd number 

dimensional 


1) dlsp"' acement = deg, ft jr m 

2) velocity = deg/sec, ft^sec or m/sec 

3) ar oration = deg/sec^, g 

4) inertial axes = deg/sec, g 
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d) generalized forces: moments and forces in y2C/v"a form 

(rotor #1 parameters, bodv axes); torque in 
form (rotor #1 parameters) ^ 
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8. UNITS 


The progran will work with English or metric (Sl) units for 
Input and output. Some of the Input parameters and most of the Internal 
program parameters are dimensionless (based on the rotor radius, the 
rotor rotational speed, and the air density), Tlie units for Input and 
output parameters are based on the consistent mass- length- time system 
(foot-slug-second or metcr-kllogram-second) , with the following 
exceptions : 

a) The aircraft gross weight Is Input In pounds or kilograms. 

b) The aircraft velocity is Input in knots for both systems 

of units (alternatively tte dimensionless speed can be input). 

c) Power is output in horsepower for both systems of units. 

The "dimensional" output for angles is in degrees; the "dimensionless" 
form for angles Is In zadlans. 
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9. AIRFOIL TABLE PREPARATION 

This section describes a program that constructs airfoil table 
files in the form required by the rotor analysis. The program will also 
print or printer-i ^t the airfoil data in the file being ci*eated or 
in an existing file. The airfoil tables are constructed using either 
analytical expressions or an airfoil table deck (in C8l format). The 
subprogram functions and namelist input labels are summarized below. 

Subprogram 

Name 

MAINTB Airfoil table preparation (main program) 

AERCT Interpolate airfoil tables 

AERCPP Printer-plot airfoil aerodynamic chairacteristics 

Namelist 

Label 

NLTABL Table and print/plot data 

NLCHAR Airfoil characteristics data 

The structure of a job to run the airfoil table prepeiration program is 
defined below. The basic structure consists of the following steps: 

1) Airfoil file definition 

2) Main program call 

3) Title card 

4) Namelist NLTABL 

5) Fcr each radial station (CPREAD f O), either 

a) Namelist NLCHAR (OPRE.AD = l) 

b) Airfoil table card deck (OPREAD = 2) 

Sample jobs are presented below. 

Create airfoil table using analytical expressions. 

DIEF FT40F001,, AIRFOIL 
CALL MAINPROG 
title card 

&NLTABL table data,NFAF»40,0PREAD=l,&END 
&NLCHAR airfoil characteristics data,&END 

?SEND 
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Create airfoil table using C8l format airfoil card deck 

DDEF FT40F001,, AIRFOIL 
CALL MAINPROG 
title card 

&NLTABL table data,NFAF=40,0PREAD=2,&END 
airfoil card deck 
?5EHD 


Print and plot airfoil table data 

DIEF FliUOFOOl,, AIRFOIL 
CALL MAINPROG 
blank card 

&NLTABL output data,NFAF=40,0PREAI^=0,<&END 
“SEND 


Ti following pages described the input variables and data for 
the airfi. il table preparation program. 


First Card 

TITLE(2( j title (80 characteirs) ; blank card for CPREAD EQ 0 


Namelist NLTABL 


angle of attack boundaries 

NAB number of boundaries , N ; maximum 20 

3L 

NA(NAB) 't. licies at boundaries, nj^ 

A(NAB) ^ at boundaries (deg, -i80° to 180°) 

Mach number boundaries 

NMB number of boundaries, N : maximiom 20 

NI‘!(N'!B) indlcies at boundaries, nj^ 

M(NMB) ’i at boundaries (O. to 1.) 

radial segments 

NRB number of segments, N^; maximum 10 

R(NRB+1) boundaries of segments (R(l)=0,, R(NRB+1)=1.) 


maximum NAB*NMB*NRB = 5000 



0PraNT(3) 


integer parauneter controlling output; EQ 0 to 
suppress; default value is i 

(l^ interpolate and print 

(2) interpolate and plot 

(3) list tables 

NMHiNT number of Mach number values for print and plot; 

maximum 10 

MPRNT(NMPRNT) Mach number values for print and plot 

NAPRNT number of angle of attack values for print; maximum 60 

APRNT(NAHINT) angle of attack values (deg) 

NFAF unit number for airfoil table file (default 40) 

OPREAD integer parameter: EQ 0 to read airfoil table and 

print data only; EQ 1 to create airfoil table using 
analytical expressions, write airfoil file, and 
print data (default); EQ 2 to create airfoil table 
using C8l format airfoil card deck, write airfoil 
file, and print data 

Neimelist NLCHAR (for each radial station; if CPREAD = l) 


CLA 

MDIV 

^LMAX 

FSTALL 

MSTALL 

GSTALL 

HSTALL 

CLF 

CMAC 

CMS 

DELO 

mi 

DEL2 

DCDDM 

MGRIT 

ACRIT 

ALFD 

GIF 


a = c, at M = 0 (per rad) (default 5*7) 

lift divergence Mach number M^^^^ (default . 75 ) 

Co at M = 0 (default 1.2) 

*max 

factor f for co (default 0.5) 

5 lUolX 

Mach number M for co (default 0.4) 
s '‘max ^ ' 

factor g, for stall c« (default 1.2) 

factor h for stall c_ (default 0,4) 

6 ^ 

cj|^ for s'-all Cjj (default 1.12) 

Cmac (default 0.) 

Cjjg (default -O.Op) 

So(defo'. 0.0084) 

(default -0.0102) 

^2 (default 0.384) 

(default O.65) 

critical Mach number at b<.= 0 (default 0.83 ) 
critical Mach number zero at e< = (default 33 •) 

drag stall angle (deg) (default 1C.) 
cd£ for stall c^ (default 2.05) 
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Airfoil Card Deck (for each radial station; if OPREAD = 2) 


I . Header 

a) Card 1, format (30A1.6I2) 

title, 30 alphanumeric characters 

NML, number of Mach number entries in Cj^ table 
NAL) number of angle of attack entries in Cj^ table 
NMD, number of Mach number entries in Cj table 
NAD, number of angle of attack entries in c^ table 

NMM, number of Mach number entries in c„, table 

NAM, number of angle of attack entries in c table 

m 

II. Lift Coefficient Table 

b) Card 2, format (7X,9F?.0); 2 or more cards if NML > 10 

Mach numbers M^ to Mj^ 

c) Card 3a, format (F7.0,9F?.0) 

angle of attack, e<^ 

lift coefficients c^^ at M = M^ to Mj^ or M^ 

Card 3b, format (7X,9F7.0); 1 or more cards if NML 10 
lift coefficients Cj^ at M = M^^p to M^^,^ 

d ) repeat card 3 for to 

III. Drag Coefficient Table 

e-g) format as or lift coefficient table 

IV. Moment Coefficient Table 

h-j) format as for lift coefficient table 

V. Parameter Limits 

a) Mj^ = 0; data extra p:< la ted for M > Mj,jj^; Mach numbers in 
sequential order 

b) = -180°, " 180°; angles of attack in sequential order 

c) NM 2, NA i 2 for lift, drag, and moment 

d) (NM+1)(NA+1) ■5' 501 for lift, 1101 for drag, 576 for moment 
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For OHIEAD •= 1, the progreun calculates representative airfoil 
characteristics using the following expressions (refer to the 
accompanying figures). 


a) Below stall 


a/^l-M^ ‘ 

M < M-, 
div 


lU, +.1 

div dav 


M -<M,, +.1 
div 




^'''ac 

'd ■ *0 

Ac^ = max (O, ) 

= max (0, (l- 


B) Stall angle 

■ '=As''=«« 


(1-M) + f,,(M-M^) 
1-M_ 




C) Stalled lift (lo<l > c<J 

o 

. r (g l»<l)co + ('“l-ot )h cfl 

Cg = signv^x) max I ^^s^^ *^s s^ s *s 

L «s**s - -<s 

max cjjf sin 2\o<\ 

Cjj = cjjj. sir 2 •< if W\ > 45° 
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D) stalled moment (|*<.\ > 

c = X sign(oc) (60- (l«x| -o<s)-75cn,f I- 1 <60° 

“ \ 60- •< 

s 

sign(<x) (90- l«l)*75cnjf + (lo^^ "60)%^. loc(>60° 

L 30 

°nif = -ic^(*<=90) = -Kc^(«<=“<^)+Cd^) 

E) Stalled drag ( lo<\ ><X.^) 

+ Cdf sin 

F) Reverse flow (lo/\>90) 

use effective angle of attack and account for moment axis shift 

o< = o< — TV sign<9< 
e 

c = C„ (z )c- + (i sincx )c, 

m m e ^ c u 


/ ^d 90 \ 

( 90 - 0 .^ y 
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